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CHAPTER 1 


INTRODUCTION 

This Is a report of the initial phase of research under Grant NSG 1128, 
awarded to the University In December 1974. The project Is part of a 
program to design a receiver for aircraft (A/C), which, as a component 
of the proposed Microwave Landing System (MLS), Is capable of optimal 
performance In the multipath environments found in air terminal areas Cl 3. 
This project focuses on the angle-tracking problem of the MLS receiver 
and deals with signal modeling, preliminary approaches to optimal design, 
suboptima I design and simulation study. This effort was integrated in 
part with the work under contract NASt-l2754-l , which was summarized in 
a final report dated May 1975 C2l!, 


CHAPTER I i . - 
SIGNAL MODELING 

The angular coordinate determination capability of the MLS has been 
tentatively ei* ^blished by a provision potentially for six (6) angle 
channels ClD* The signals in all channels have similar forms, hov/ever, 
and without toss of generality, this study has employed the azimuth 
channel as a focus to provide concrete Interpretations of the models and 
results developed. Parallel results for other channels v/ould obtain 
from suitable adjustments of parameter values. 

The coordinate system in use is the right-handed system shown in 
Figure 1-1 with the origin located at the center of the azimuth transmit- 
ting antenna. The relevant variables r, 0, ^ are re ge (n.mi.), azimuth 
(degrees, positive clockwise looking down) and elevation (degrees, 
positive up). Let p(») and respectively the antenna selectivity 

and scanning functions, be defined as follows: 

p(0g) = the far-field strength of the transmitted signal (field) 

relative to the bores ight signal strength at the same range, 
where 0^ is the azimuth angle of the observer relative to 
the antenna boresight axis, and thus p(0) = I. (2-1) 

0^(t) - the angular orientation in azimuth of the antenna boresight 
axis as observed at time t at the transmitting antenna 
site. (2-2) 

The definition of p(*) above does not recognize explicitly a dependency 
on observer elevation, which Is present due to phased array radiation 
pattern properties. The effect is recognized implicitly however, with 
notation suppressed unti 1 cross-coupling of the various angular coordinate 
estimators (and the DME) is considered. 

Further, letting 

t, = the time (as observed at the antenna site) that the 
k 

kth TO-FRO antenna scan begins; (2-3) 

T - the duration of the TO-FRO scan as observed by a 

stationary observer, (2-4) 
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Figure ll-l Coordinate System 
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then is specifically defined by fhe antenna scanning motion as 

a function of Its argument on the Interval Ctj^^tj^+Tlc 

The coordinates of the receiver (the A/C) at time t are r(t), 6<t), 
<{)Ct), A di rect-path-slgna I event arriving at the A/C at time t Is one 

P { "t* ) 

that necessarily left the transmitting antenna at time t — - — , where 

c is the propagation velocity of light. Hence, the component of the 

received signal observeo at time t and due to transmitter emissions 

P ( "t” ) 

derives from transmission occurring when t, < t — ^ < t, t T. This 

component Itself comprises two contributions Yq^T)# follows, 

due respectively to direct-path and reflection-path propagation: 

YpCt) - a(t)pCe^(t - ) - 0(t)3 cosCdi^(t - ) t (2-5) 




( 2 - 6 ) 


where 


and 


Cj (t) rj Ct) 

(t) - 0^(t)pCe^(t ^ cosCw^Ct ^ ) + Pp ] 

(2-7) 


a(t), = amplitude and phase of the direct path component (2-8) 

(ij = radian frequency of propagating r-f energy (2-9) 

0 

r^ (t),0p (t) = range thru and azimuth of the 1th reflection point (2-10) 

1 j 


ai(t),pR 


= amplitude and phase of the component reflected by 
' the ith reflection point (where, for each I, a. (t) = 0, 

Iff ^ ^ 

Under the assumption of isotropic reflection at the Ith ref lection point, 
located at r-r (t), 0 r. (t), <{>R,(t), the range through reflection, Vj (t), 
is easily shown (see Figure 11-2) to be 


Ti 


(t) = rRj(t) t rflD.(t) 


ARi 


(2-12) 
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where 


^AR 


(t) = Cr^(t) - 2rCt)rD + rD^(t)H 


:/2 


M 




(2-13) 


and 


?(t) = cos^iCt) cos^ip (t) cosCeCt) - 0p (t)H + sin4>(t) sln4ijj (t) (2-14) 


R 




In view of the generally independent motions of the A/C and a reflector, 
the equations (2-7) and (2-12) through (2-14) are clearly approximations. 
A more exact expression for, say t , the time of transmission, given 

/S 

receipt at time t (I.e., t^ would be the argument of 0^( ) In y^ (t) 
generally) would have the form ^ 


^AR, ^R. '^AR. 

t = t - -J L (t - — ) 

c c c * 


(2- 1 4a) 


based on an argument similar to that preceeding (2-5) above. On the 
basis of the linear approximation 


AR. 


'■r!+ - 


i(t) 


AR 


) = (t) - (t) 

t^i R. 


i(t) 


the expression for t becomes 


''ar''*'’ '■r.'+>'‘ar /+> 


(t) rp (t)r (t) 

= t ! + _ ' • ' 




The original approximation for t (in <2-7)) results when the term 

A 

Td ’“flo /c^ is dropped, which seems reasonable under the circumstances. 

Kj nKj 


In the strict sense the expression in (2- 1 4a) does not give the exact 
value for t , but is adequate to show the approximate magnitude of error 
in the original approximation (in (2-7)). 
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The basic model has now been presented in terms of a number of 
parameters potentially time-varying within a scan Interval In an un- 
specified manner. The update rate of the MLS (13.3 Hz in Azimuth) Is 
sufficiently high relative to the dynamics of the environment to allow 
simple approximations for these parametric functions within a scan 

interval. For < t < t. , (which contains the kth scan interval, 
rCt) ^ ~ 

tk < t — < tk + T, at the A/C), the approximations used are the 

following, which seem reasonable: 


a(t) = aCtk) 

e(t) = 0(tk) 

r(t) ^ r(tk) + r{tk)nt - 




(2-15) 

(2-16) 

(2-17) 

(2-18) 

(2-19) 

( 2 - 20 ) 


it is emphasized these approximations are limited to a given scan. Clearly, 
the result in (2-12) through (2-14) generally Invalidates (2-20) In any 


but a local sense, with rj^(tk) determined even then by (2-12) through 
(2-14), (and r-r (t, ) by the derivative of the expression In (2-12), as 

' I K 

was done in the simulation study; see Section IV). 


Under the assumptions (2-15) through (2-20) the received signal 

can be expressed in terms of a time variable local to the scan being 

r(t) 

received at the A/C. First, note that t ^ — , the argument of 

in the received direct path component, Vq(+>, (2-5), as a direct result 
of (2-17), can be written as 


r(t) - I _ rC^k) j _ rC^k) ^ ^ rC'^k) 
c e c k c 


( 2 - 21 ) 


Consequently, the "received scan interval," defined heretofore by 
r(t) 

t, < t - ■ ■ ■■ I < t. + T, can be described as fol lows; 

K C *■ K 
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( 2 - 22 ) 


r(t. )/c r(t. )/c j 

k I - rCt|^)/c - - k 1 ~ r(tj^)/c I - r(t|^)/c 

and thus on th© kth scan Interval we may write In general 


rCt. )/c 


(2-23) 


where t Is a local time variable on the kth received scan interval, vary- 
ing as fol lows: 


0 < T 


- 1' - Ktj^)/c 


(2-24) 


with the same scaling as t, I.e., dt/dx = I, (note a change in the 
duration of the received scan period In (2-22) and (2-24) due to A/C 
motion). The quantity 


. A. , 

\ I - 


(2-25) 


Is the time of reception at the A/C of a signal event transmitted at 
time tj^, in particular the start of the kth scan. Errors in determining 
the arrival of this event at the A/C translate Into errors In starting 
1he local x-clock. It is th© function of the Barker code, however, to 
keep these errors small. 

As a consequence of the above development 


t 


r(t) 

c 


= t, + (I 
k 


r(t^) 




(2-26) 


and, using approximations (2-15) and (2-I6), the direct path component 
y^ can be redefined as a function of a global discrete-time variable 
t|^ and a local continuous-time variable x, as follows: 

yjj(tj^,x) = a(t|^) pCe^(t|^ + x) - e(t|^)Ii cosCwQCtj^)x t e(t|^)] (2-27) 
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in which the -r(tj^)T/c term of t2~26) has been dropped from the argument 
of because of Its relative Insignificance, and w^tt^), S(tj^> are 

defined as fol lows; 

* j^tD (I - r(tj^)/c), as the apparent r-f carrier frequency 
WoCt.) = ( . ‘ (2-28a) 

l^w.p - m rft. )/c, as the apparent t-f carrier frequency 
‘‘‘ ^ (2- 28b) 

0(tj^) “ Bq + initial phase on the kth scan period (2~29) 


The expected doppler shift (neglecting relativity) is present In (2-28). 

The phase parameter 6p In (2-5) and (2-8) was modeled originally 
as an independent random variable uniformly distributed on C-Tr,irJ, 

Equation (2-29) above for suggests that 6(t|^) then is also a 

random variable but one possible highly correlated with 0p. As a prac- 
tical matter, however, the Increment ^ large multiple 

of 2 it radians and can nevo-r be determined accurately enough to make 
valid use of the correlation implied by (2-29) above, for two reasons: 

1. The t|^ are determined from the Barter decoder output via (2-25), 
and with a noisy Input the uncertainty may become a significant 
part of a bit pulse width (approximately 67 ps, also many periods 
of the 5000 MHs r-f carrier). 

2. The quantity t^^ - tj^^l is not constant from scan-to-scan, and 
may even be jittered to suppress propeiier modulation effects. 

On this basis the B(t|^) also were modeled as Independent random variables 
uniformly distributed on C”'r,'r3, which seems reasonable. 

The ith reflection component, y_ (t), (2-7), can also be expressed 

"i 

in terms of the local receiver time variable t, associated above with 

the direct path signal, Under the approximation of r-p (t), (2-20), we 

I 

may write 


t - 


r^ (t) 
’i 
c 


t(l 




) - 


c 


+ t. 


r(tj^) 

c 


(2-30) 
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Substituting from (2“33) for t and (2-25) for t^^ gives 

Ty (t) r-r (tj5 ) ^ 

' t ' t '^k * 

t = t + -cCI - Ct.)/c) 

C C 1 ^ 


(2-51) 


Ar, (t, ) Ar, (t, ) r(t, )/c 




r 'k 


r ‘k' 


where, in generj! 


-pr-[y7j+ •'(I - rx'V/c) 

(2-32) 


Ar. (t) = ry Ct) - r(t) 
' 'l 


Ar. (t) - Py (t) - r(t) 


(2-33) 

(2-34) 


Scan data Is received through the reflection path, clearly, when 
I , , . '■'V'= ,n - 

^ - I - fj(+|^)/C c ‘ c ’ ' I - ■ 


c 

(2-34a) 


Neither of the terms involving the time derivatives Apj Ctj^) and r^ (+[^) 
in (2-32) is important In the argument of hence, using the approx- 

imations (2-18, 2-19) for a,(t), 6„ (t) respectively, the ith reflection 

component y,, can be redefined with bivariate argument (t.,T) as follows, 

Ki T ^ 

0 < T < J - = 


I 




^rj(t-) 

a,(tj^) p :e^(t^ -4- T - ^ - 0 r^CV3 


Ar.(t. ) ^ 

“='^“0, '+k>^' -V— i i 

* K 


0, 0 < T < Ar. (tj^)/c 


(2- 35a) 
(2-35b) 


where 


“o.'V = 

I 


^(1 “ r j (tj^)/c), as r-f carrier frequency (2-36a) 




(tj^)/c, as i-f carrier frequency (2-36b) 
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; r 
i I 
1 


«t‘\’ ' ®R, + “cC+R^ - 

= ®R, ^ “A - : 


(2-57) 


(2-37a) 


Ar.(t,) Ar.Ct ) rCt )/c 

= \ * “cl^+k - • T ■ - <2-37b) 

I k 


A more compact form of (2-35), which will be useful, employs a bivariate 
argument in a., as follov/s; 

Ar. (t, ) 

Y„^(1-|^,T) = + T - —1^ ) - 6^^(t|^)3 


where 






“I'V' Ar,(tk>/c <T< , _ 


0, 0 < T < Ar j 


(2-38) 

(2-39a) 

(2~39b) 


i / 


It was established, following equation (2-29) for the 3(t|^) in the 
direct path component, that the 3(tj^) were independent random variables. 
The 3(t|^), however, are functions of properties of the transmitted 
signal which the 3.(tj^) In the reflection components share. For each I, 
the initial value of 8j(tj^), say pj(tj), was modeled as an independent 
random variable uniformly distributed on [-ir,ull, but as a consequence 
of the above, it was felt that the scan-to-scan change in 8j, e.g., 
Si(t|^^|) - 3jCtj^), should be related to the corresponding change in 3, 
i.e. 3(t|^_l^j) - (2-29) and (2-37) gave the relations 

CD 

- 6(+^)3 - ^ Ur (t„ ) - Ar (tR )1 

k+i k 

(2-40) 


(j 


- C3Ct,,^r) - 3(t, )J - ~ Ct,.., - t,)Ar,Ct. ) 


k+I 


k+1 


k' 

( 2-40a ) 


in the latter of which relatively insignificant terms corresponding to 


t 


the last in (2-37b) were dropped. In the simulation this relation assured 
the smooth evolution of the simulated multipath interference phenomena. 

Specializing more concretely now to the signal In the (linear) l-f 
channel of the receiver, where the total signal, y, comprises direct 
path, reflection path and receiver noise components, we assume the i-f 
bandpass characteristic is appropriately symmetrical about the nominal 
bandcenter frequency, The noise then, a process generated at the 

receiver and hence modeled with sample function pCt), Is assumed to 
be a stationary, bandpass, zero-mean Gaussian process with variance 

ft quadrature expansion of the noise n(x) is possible, therefore, 
and the following formulation of the composite signal y with bivariate 
argument results: 

y(tj^,T) = y^Ctj^,T) + yp^(tj.,T) + n(x) (2-41) 

where 

yj^(tj^,x) = a(tj^)pll0^(tj^+ t) - e(tj^)]cosCwjpT+ ($(tj^) 

" '^D f\A)sintO|pT 

c. ^ 

Ar, (t. ) 

= p, (+k,T)pC6^(t, T - 

COsCtOjpT + (e.Ctj^) - (tj^)/c) 3 

I 

— yj^ (t|^, T )cos[i} I pT — (t|^ , T ) s i nu) I pT 

c s 

n(x) = n^(T)cos(D| pT - ng(T)sinW|pT (2-46) 

in which 

“ 9Ctj^)-lcosCe(t^) - w^Tr(tj^)/c] (2-47) 


(2-44) 

(2-45) 


- oj xr(t. )/c)3 
C K 

(2-42) 

(2-43) 
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Yp " 9t+|^>3sinCp(t|^) - a)^Tr(tj^)/cl 

(2-48) 

Ar, (t, ) 

vr <+,<.•') = I“i<V'’pi^V+k ^ ■' — s — ’ - V‘+k’^ 

c I i 

cosllpjCtj^) - (t|^)/cl] (2-49) 

Ar.(t.) 

Yr^C^t) = Ia.(tj^,T)p:e^(t^ + T ^ ) - 0 R.Ct,^): 

sinC3.(t|^) - w TrT (t, )/c3 (2-50) 

IK C I I K 

and 

n (t), n (t), the quadrature components of the noise, are independent, 
0 s 

stationary low-pass Gaussian processes, each with mean zero, variance 
0 p^, and two-sided bandwidth equal to the i-f bandwidth. (2-51) 

Hence 

y(tj^,x) = m^Ctj^,T)cosw j j-T - mg(t|^,T )sinwj pX (2-52) 

= M(tj^,x)cosCw j pT - r(t|^,x)l (2-53) 

where 



= Vp + Yr %.r) + n (T) 

c c 

(2-54) 


' Vp (+|,,t) + Yr + n^(T) 

s s 

(2-55) 

M(t, ,x) 
k 

= '^mc^(tj^,x) + 

(2-56) 


= arc tanCmg(t|^,T)/m^(t|^,T)3 

(2-57) 


The principal signal modeling results to date have been presented. 
These results have been used in -rhe digital simulation studies described 
in Chapter IV and in the optimal receiver design analysis of Chapter HI. 
in the simulation the local time variable x was discretized, of course, 
a sampling rate equal to the I-f bandwidth (160 kHz) being chosen. 


13 


other small refinements, such as dropping certain r/c terms, etc., 
were made also where permissable; reference Is made to Chapter IV 
for a more detailed description of the simulation study. 

The optimal design studies undertaken to date and reported in 
Chapter III were specialized to a signal model Involving only the 
direct-path and receiver noise components. This Initial focus on 
the mu Itipath-f ree case was justified on several counts: 

1. The study establishes for reference purposes the optimum 
level of performance obtainable 'n a multi path- free 
environment. 

2. It provides Insight into the structural properties of MLS 
receivers resulting from the requirement for optimal 
performance in a very basic and generally applicable 
disruptive environment. 

3. The algorithms that result under this specialization represent 
lower bounds In computational complexity in relation to those 
optima! processors associated with the more complex environ- 
ment of Interest. 

The first point cannot be overemphasized; the MLS receiver will operate 
most of the time without significant multipath interference, and its 
performance during such periods Is not outside our concern in this 
study. 

In re-extendIng our concern to the multipath-corrupted signal, 
the results obtained above will probably not be used directly because 
of the very high dimension of the parameter space associated with the 
total reflected signal 

where y^ is given in equation (2-38), Instead a stocastic process 

model (really a random field model) with associated sample function 

yo(t, ,t) is being considered. This is done in conjuntlon with considerations 

of the estimation algorithm to be used, since the statistics of interest 

14 ■ 


depend to some extent on the latter. Second-order statistics on the 
scan interval and local to the arrival time of the direct-path pulse 
are expected to be Important; variations of these statistics with 
are the basis for an adaptive approach to the suppression of multipath 
interference. If first- and second-order statistics turn out to be 
sufficient, a unique Gaussian random field model is implied and the 
possibility of an adaptive algorithm of manageable dimensions. 

Lastly, a state-var table model Is needed as a prerequisite to the 
application of modern recursive estimation techniques to this problem. 
Such a model generally would describe the (short-term) evolution of 
the A/C angular coordinate and other parameters of interest and relate 
these to the observed receiver signal. In the desired model the "state" 
should Include the A/C angular coordinate, of course, but also all un- 
knov/n parameters whose estimates are needed for use in the angular 
coordinate estimate calculation. In the locally optimum estimation 
scheme to be described in Chapter 111, for example, the latter parameter 
set included the signal strength parameter a(t) and the apparent i-f 
carrlfT frequency ^^(t) (Note: the latter estimate may not be needed 

if the nominal i-f carrier frequency is sufficiently great with respect 
to the doppler shift (which is about 800 Hz/ 1 00 kts speed, for a 5000 MH 
r-f carrier frequency)). In addition, in an adaptive design an "aug- 
mented state" may appear and include some parameters associated with evo 
lution of the "state-proper," thus allowing some estimation of the 
mode! as well as the ^variab I es of prime interest. The forms of models 
needed In the recursive estimation study will be determined later when 
recursive estimation algorithms are considered in greater detail. 
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CHAPTER 1 I I 


OPTIMAL MLS RECEIVER DESIGN 

Under the optimal tracking criterion two approaches to receiver 
design were considered: 

1. Locally optimum estimation C3D in which it is assumea the 
tracking error Is always small and constant on the interval 
of observation (the scan Interval), and thus can be esti- 
mated as a parameter without requiring a knowledge of the 

a priori statistics of the A/C angular coordinate of interest. 

2. Recursive state estimation [_ 4 , Chapters 7, 8, and 9] In which 
a state space model is assumed and used to provide smoothing 
of the evolving estimates and also extrapolation into and 
possible through periods of signal fade. 

Quadrature detection of the l-F signal arises in these algorithns, and 
to preserve the option for analog implementation of these stages con- 
tinuous-time formulations were used mostly. The principal signal 
model used in the Initial studies was the direct-path-component-p lys- 
noise specialization. 

Locally Optimum Estimation 

In this formulation the received signal on the present scan interval 
is modeled as an observations process {y(t),t e C0,T]> with sample 


function of the general form 

yCt) *''*s(t,m+e,a) + n(t), t ^ Co,Tl] (3-1) 

where 

t = local receiver time on the present scan Interval (3-2) 

C0,Tj=the present scan Interval wrt local receiver time (3-3) 

s(t,*,»)= a sample function of a process or a sure function with which 
the modulation or parameter of interest is associated. (3-4) 
n(t) = a sample function of a white Gaussian noise process (n(t), 

t E C0,Tll} with known power density N^. (3-5) 


m+e = true present value of the unknown parameter (or parameter vector) 
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•to be estimated, assumed constant for t e 


(3-6) 


m = estimate of the unknown paramster, based on observations 
through the last, but not the present, scan Interval 
(possible extrapolated to the present, however) (3-7) 

e = error in the last estimate, m (3-8) 

a = a random parameter (or parameter vector) in which we have 

no estimaticn interest (3-9) 

To avoid technical difficulty with mathematical processing of white 
noise the integrated observations process {Y(t),t e LO,!]} is con- 
sidered instead of {y(t), t e C0,T3}, where 

Y(t) = /Js(x,m+e,a)dT + N(t) (3-10) 

and N(t) is a Wiener process for which 

^N(t)N(t)^ = Nj^min(t,T) (3-11) 

(and where ^ denotes mathematical expectation). 

The objective is to estimate the (unobservable) error e and to 
refine the parameter estimate m using the error estimate. Following 
Murphy C3H, probability measures relevant to this objective In a 


locally optimum estimation context are the following: 

Q = measure corresponding to {N(t), t e C0,T^} (3-12) 

S = measure corresponding to {s(t,mte,a), t e C0,T3) (3-13) 

Y = measure corresponding to {Y(t), t s C0,T3 |s} (3-14) 

s 

P = measure corresponding to {Y(t), t e C0,T3} (3-15) 

A ■ 


An estimate e of e is said to be locally optimum C33 a+ e == ©q If two 
conditions, stated below, are satisfied; all integrations are taken over 
the underlying elementary event space (universe), denoted fl: 

I.) The estimate e is locally unbiased at e^, I.e., 

a) /q<© - ~ ^ vector), and (3-16) 


(3-17} 


b) I— f^(e - e^j'^dP- = Kthe unity matrix), 

3e 0 ®0 

T 

where ( ) denotes the transpose of ( ) . 

2.) The estimate e has minimum mean square error among all locally 
unbiased estimates of e at e^, l.e. If a Is any estimate of e 
locally unbiased at e^, then 

/(e - e^lCe - eo>''^dPeQ < f(o - 

In the sense of the usual "weak” order relation on a 
set of non-negatIve definite square matrices of the 
same dimension. 

A 

Murphy C33 has shown that the locally optimum estimate e of e at 
e = 0 is given by 

e = ^q'Aq(Y) C3-I9) 


where is the vector whose Ith component, A^ , is given by 


A^m =' 

I 


dP 

( 2. ) 

3e. n dQ 


dP 




e=0 


, otherwise 


(3-20) 


and 


*0 = 


(5-21) 


in which Yq is given by (3-10) with e = 0. if e = 0, the residual mean 
square error (l.e., the error covariance matrix) associated with this 
estimate is that Is 



Radon-N I kodym derivatives associated with this estimate are the fo I lowing 


dY 

c 

w 


- exp ( ■jli' /Qs(t,m+e,a)dYCt) - /QS^{t,rnte,a)dt) (3-23) 


I rT„2. 
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(the first Integral above being a Wiener integral) and 


dP dY 

„ - /’■g^ S(ds), 


(3-24) 


this latter integral being taken over the space of realizations of 

{s(t,m+e,a) , t e CO^T^}. When s(t.,m+e,a) - s(t,m+e)/ t.e. is deterministic, 

then S Is degenerate, I.e., there exists a realization of the process 
“ dPe dY„ 

{s(t,m4*e), t e H0,T!]} such that S({s>) = I, In which case . 

Another special case of interest Is when s(t,m+e,a) is a deterministic 

function of random parameter a, which ranges on the reals, R|, with 

induced probability measure F. Then the collection of real izations of 

the process {s{t,m+e,a), t e C0,TD} is indexed by a e Rj,and we may 

write 


dP dY (t,m+e,a) 

dr=/R,-V— 


(3-25) 


where dYg(t,m+e,a)/dQ denotes the form in (3-23) with the functional 
form of s(t,mte,a) substituted as required. If y is Lebesque measure on 
and P is absolutely continuous with respect to y, then we may also 
write 


d P dY ( t , m+e , a ) 

e _ c s' ' 

d^ " JRj dQ 


dP 

-r- (a)dy(a) 
dy 


(5-26) 


where 


dp 


(a) 


is the probability density of the random variable a. 


Application of locally optimum estimation to the MLS receiver 
design problem, specialized to the direct-path-signal-plus-noise case, 
is made by letting s = y^^, i.e., referring to (2-27), 


s(t,m+e,a). = apCs^i^Ct) - eljcosiluQt t $3, t e L0,T3 
and making the associations 

V 

m+e - I a 

\ 



(3-27) 


(3-28) 


9 


a = B 


C3-29) 


It is illuminating to regard 3 as a known parameter temporarily and 
regard s( ) as a "sure" function of t and the parameters e, a, oiq. 

In this case, as explained following (3«24), the measure S is degenerate 
and 


dP 

e 

dQ 


dY C ) 
s 

"3(5 


where 


expC ^ rpCe,(t) 


- eI]cos[;wQt+ BlIdYCt) - 

(3-30) 


q|(8) = * elcos^Cw^t + e3dt 

' /Jp^Ce^ct) - eUdt 


(3-31) 

(3-32) 


Then, from (3-20) 


/ “«/JpC0/\( + ) - 93cos[wQt + e3dY(t) - |^qj(0) \ 




/JpC0Aft) " 03cosCw(jt + s3dY(t) - J q|(6) 




(3-33) 


QpC0A<t) - e3ts Induct + B3dY(t) 


y 


where 


A A <\ 


0, a, [Oq are last estimates, corresponding to m 


p(e^(t, - e) 


0^ = (0 Jt) - 0) 
e A 


. - A 

qi(0) = — “ 
de 


-z/JpCeAtt) - 03pC0A(t) - .03dt. 


(3-34) 

(3-55) 

(3-36) 


The last quantity, qi(6J, is very nearly zero (hence qi(9) is very nearly 
a constant), except when 0 is near the ends of the Interval of coverage, 
e,g. (-60®, +60°) for azimuth. To an extent the bivariate nature of 
p,( ), i.e. the channel -cross-coup I ing effect discussed following (2-2), 
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also affects the functions qi( ) and qi( ). Under the assumption these 
combined effects are negligible, the approximations 


qi(0) = constant (=/QP^C6y^(t)ddt) 

q^Ce) = 0 

are made with the result that 


where 


M = 


/ A , A A 

~apCe^(t) - eDcosCcopt + 3I 

A A 

pCB^ct) - eUcosLwQt + 0] 




tpCe^ct) - eDsinCw^t + ell 



(3-37) 

(3-38) 


(3-39) 


(3-40) 


The result (3-39) follows from setting e = 0 In Y(t) In (3-10) and then 
substituting the resulting Yq Into (3-33) and simplifying. 

The above implies, then, by (3~2| ) that 


/ C/QMdN(t)X/QM‘dN(t)]dpQ 
0 g 

= <([/jMdN(t)X/jM'''dN(t)]^ 


(3-41) 

(3-42) 


which, since {N(t),t e C0,Tl]} is a Wiener process with ^NCt.hKr)^ as 
given In (3-11), can be written 



(3-43) 


Substitution from (5-40) and use of some trigonometric identities (and 
some nearly obvious notatlonal shorthand) gives 


' Ng ^O! 


t^pZcli-SSiLI) ia+pp(SiJ^, 


\ 


-; pf ( i ±. p 2 u 


) p 2 ( J _ L | 252 U , - i + p2 , lin | U ) 


d+ 


btpp( gJ j f . J -) -;tp2(iiB|0) ;2t2p2(i-:-S2£2LI)j 


And asymptotically as uq 9®'*‘s very large), 


(3-44) 


/ «"q2 
^ ~ ir ~ 


V 


0 


a_ 

2 


A A 

a^qsCe) 
— 2 


C3-45) 


v/here, under the approximating assumptions made. 


q|(6) = /Jp^C0^{t) - e3dt 




q, = 0 


T 

= /QP^E0^(t)3dt = q., a constant 


C3-46) 


= /!p^L9fl<t)ldt = q,, a constant 
q, - 0 ^ ^ ^ C3-47) 


q^O) = / p^Ce^Ct) - eHt^^dt 


q, ™ 0 




(3-48) 


.-I 


Consequently, the covariance matrix, , for the estimation erroi , 
e - 6 q, Is approximately 


* ~ N 

'0 0 


a/Ctt^q^) 
0 
0 


0 

2/q| 

0 


0 

. 0 


(3-49) 
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and the locally optimum estimate e is obtained by combining (3-46) 
through (3-49) and C3-33) with (3-19): 

■ »3cosC<«ot + B]dY(t) \ / ®e \ 


e = 


/JpC0;\(t) - e^cosCiUQt + e]dY(t) - a 


or, alternatively (recall g was assumed known) 


/ 


/" " 03coswQtdY(t)) 

- sJna(/J!'pCe. (t) - sDsintij„tdYCt))I] 


e = 


-a+ — Ccose(/^pC«y^(t) - e!]coswQtdY(t) ) 

- sln3(/QpC9^(t) - eJsInujQtdYCt) )H 
^‘^jCcos3(/JpC0^(t) - 0:}tsInwQtdY(t)) 


6lq,(0, 


"uO 


f 

/ 


+ sInB(/]|pC0*( + ) - eljtcosto^tdY(t))] 


/ 


(3-50) 


(3-51 ) 


For Instrumentation purposes, dY(t) would be replaced with yCt)dt and 
physical dumping Integrators employed to integrate the six outputs of three 
v;elghted or gated quadrature detectors, corresponding to the above. Other 
significant observations that can be made Include: 

1. The estimator Is dependent upon the prior estimate values 
but not the noise variance parameter, 

2. The error covariance Is proportional to noise parameter Nq, and 
of particular Interest Is the variance associated with e„, the 
estimated A/C angular parameter: 

<®e> ^ ^ > <3-52) 

where the first factor is a measure of the ratio of noise-to- 
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estimated signal power, and the second factor Is the reciprocal 
of the Integral of p^, given in (3-47). 

A 

3. The estimate e In (3-50) involves an Integral less the 

C( A 

prior estimate of a, a. Hence the integral is a global esti- 
mate of oi, given y(t) and prior estimates of 0 and The 
error covariance of the estimate Is proportional to but 
independent of other parameters (under the stated assumption 
q|(0) ~ 0 for all 0 of Interest. 

Returning now to the more realistic situation in which p is regarded 
as a random variable uniformly distributed on the interval C-HjirD, as 
described following (2-29), we have from (3-30) 

dY (t,m+e,a) 2 t 

= exp(- — )exp( ^ -e^cosCojQt t pHdYIt)) 

= C.expC (i^ (Y)cos3 - U (Y)sin3l| 

e IN-. Cg 

ig (Y) 

= CgexpC ^ t i|g(Y) cos(B - arc tan ( ^ ^^yy )>l 


where 


C ^ exD ( - ) 

^e ^ 4 Nq ' 


(Y) = /JpC0^(t) - ejcos ftiQtdY(t) 
e 

l^e^Y) = /JpCe^(t) - e]sino)QtdY(t) 


(3-56) 

(3-57) 

(3-58) 


and pjO) is as given in (3-32). Hence from (3-26) 


I 


dP C 

dO^ ' 27 I •'1^ + l| - arc tan( 

” —TT Or e 


S 

6 

'c 


(3-59) 


= C \A ^ /[2 (Y) + 1| (Y)) = CJA ^ E^(Y)) 
e 0 No Cg Sg e 0 e 


(3-60) 
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where 


E^(Y) ^ ^\'i CY) + \'i a) (3-62 

e c s 

e e 

and Iq ?s the modified Bessel function of the first kind, zeroth order. 
The above relates to other results, which are well known C5D.C63. 

Substituting the above results In C3-20) we obtain, after some 
manipulation, the following: 


; , ''o ,T 

4N„ - N„ 1q^(Y) • EqCY) Jo 


i.fV) T-; . 

^ /oPcoscgtdy 


I, (Y) I (Y) 


^ 0 0 rT* . * 4 ..VN 

L vn * E„(Y) 


Aq(Y) = 


aq , C 0 ) 


1 . (Y) i 


^ 1 . 1 / 0 0 fi: 

■2N^ Nq ^ Iq (Y) * EqCY) 


» Q Q ^ 

I (Y) ■ E (Y)' 

Ug U 

1, (Y) I (Y) 

a / 0 0 fl''. . " , .w 

- c ^ iriYr * /op+"‘^V^y 


i. (Y) 1 (Y) 


U (Y) * -E°4¥-) 


In which new notation appearing Is defined as follows: 

/V 

P = pCs^ct) - e] 
p - pC0^(t) - 0“1 

'o - 'o‘|r 


E^(Y) = /If 
0 c^ 


per (3-62) 


1^ (Y) = J^pCo^It) - eHcoswQtdY(t) per C3-57) 
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J 


Ig tY) = /JpCe^(t) - eDsinw^tciYCt) per (3-58) 


(3-69) 


, (Y) ^ 1 ^ E„(Y}) 

'o ' “^0 ° 


(3-70) 


where Ij is the modified Bessel function of the first kind, first order. 
A simpl If Iciiii'Ion of (3-63) Is possible with the definition of four more 
quantities similar to (3-68) and (3-69) above; 


(Y) = /JpCe^(+> “ 6>oswQtdY(t) 

Jg (Y) = /^pie^Ct) - eHsInioQtdYit) 

(Y) = /JpCe^(t) - e^tcosw^tdYCt) 


(3-71) 


(3-72) 


(3-73) 


K (Y) = /IpC6Att> - e3tsInu„tdY(t) 
Sq u a u 


(3-74) 


Further, the quantity I, (Y)/U (Y), essentially a soft-lImlter function 

0 ^0 

with Initial slope = 1/2 and maximum values of + I, will be denoted by 
Lq(Y); it Is closely approximated by an arc tangent expression, I.e., 

A 

L (Y) t = I arc tan C i . E„(Y)3 (3-75) 

Nq 0 ^ 


with these definitions Aq(Y) can be written 
/ .tt^q,(e) aL-(Y) 

/ - -1C COYT ^'c ^ 's <'^>3 


Agm = 


aq,(6) L (Y) „ 

- -4 + mV r V V C r (Y) + i: (Y)] 




aL (Y) 

jrrjyfD^ CY)K (Y) - 1^ (Y)K (Y)Ii 

Vo 0 ®0 ^0 '^o 


(3-76) 
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which together with C3-67) through (3-69) and (3-71) through (3-74) 
indicates the amount of processing of the observation y(t) (or Y(t)) that 
is required by the optimal processor which Is denied the value of the 
carrier phase parameter, 3. Further, a comparison of (3-33), the 
known 3 case, with (3-63) apove shows the two expressions for A are 
analogous If the following interpretations are made: 


cos 3 ” '■ [” 


(Y) 


(Y) 


IcqCY) 


= Lq(Y)( 


CO 


(Y) 


Eq(Y) 


(3-77) 


sin 8 = 


i, (Y) u (Y) 




(Y) 


1 fY)f 


“0 


(3-78) 


These may not be optimum or even good estimates of cos 3 and sin 3 in 
the usual sense, but locally optimum estimates of 0, a, oJq result from 
their use. Certainly it can be shown that for very large signal-to- 
noise ratios 


cos 3 ■* cos 3 


(5-79) 


stn 3 


sin 3 • 


(3-80) 


In general we settle for less, however, noting, for example 

(cof 3)^ + (sin^3)^ = Lq(Y) (3-8!) 

which is less than unity for finite signal -to- noise ratios. 

This is as far as the study of the locally optimum estimation 
technique has progressed. The computation of the error covariance 
matrix for the random 3 case has yet to be done (if further effort 
on this direct-path-slgnai -plus-noise case is warranted). This matrix 
is not a function of the observation Y(t)j it is needed to complete 

A 

the solution for the error estimate e via (3-19) as well as providing 


the desired reference performance measure. The structure of a representa- 
tive Implementation of this receiver algorithm Is shown in Figure IM-I. 

The quadrature detection represented by the calculations of 1^ , 1^ , 

Jsq# Ksq, (3-68), (3-69) and (3-71) through (3-74) respectively, 

would be done with analog circuitry, and the remaining calculations 
associated with the error estimation, (3-i9), (3-76), (3-67), (3-75), 
etc., as well as the A/C angle estimate update, would be done in a 
digital microprocessor. This algorithm requires a knowledge of the 
noise power density Nq. In a baseband simulation of this algorithm the 
quadrature detection integrations, (3-68), (3-69) and (3-7!) through 
(3-74), would be represented by integrations of the low-frequency 
(i.e. difference frequency) components of the Integrands shown In the 
equations. 

The extension of the locally optimum estimation algorithm to 
multi path-corrupted signals will be an early objective In the continuation 
of the project. The problem wilt center on modeling as a process the 
reflection component Yp(T) of the received signal y(t) = yp(+) + 
described In Chapter !!, and then calculating the required likelihood 
ratio (Radon-Ni kodym derivative). Kai lath’s work C63, H7], will 
help in this second aspect. The successful extension to the multipath- 
corrupted signal is expected to result in an adaptive algorithm, which 
approach is discussed further in the conclusion to this chapter. 

The complexity of the optimal quadrature detector processors that 
liave resulted motivates a serious consideration of the rectifier-type 
of envelope detector and the processing optimally of its output to 
produce the estimate of the A/C angular coordinate. The difficulty 
in considering this option from the locally optimum estimation view- 
point stems from the possible inapplicability of the latter model (as 
formulated by Murphy) fo the recti fled'^enve I ope process and the possibility 
that a suitable likelihood ratio expression cannot be found in the 
literature and will have to be derived. These problems are related to 
the fact that, even in the multi path-free case, the corruption of the 
envelope process (induced by receiver noise) has the following properties; 
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Figure 1 1 t-l 


Structure of a Locally Optimum Estimator Receiver 
for the Direct-Path Signal + Noise Case 















1) The corruption Is not additive, 

2) The corruption Is not Gaussian, 

3) The corruption does not have mean value zero. 

The rect I fler-enve lope-detector does seem to offer some economic advan- 
tages In Implementation, however, and the associated analysis and 
design problems will be studied. Concern to date for these types of 
envelope detectors has been somewhat informal and is discussed in 
Chapter IV on suboptima I design. 



Recursive State Estimation 


Recursive state estimation (RSE) has some appeal for the MLS 
application, but algorithms of this class In general have some undesir- 
able characteristics which must be understood and whose effects minimized. 
In addition, general RSE algorithms are based upon system models In 
standard state variable forms, and algorithms considered for the MLS 
application must accomodate the special form and uncertainties in the 
MLS signal model. The adaptation and evaluation of recursive state 
estimation techniques for the MLS application Is presently at an early 
stage. This is only a brief description of some relevant properties of 
recursive state estimators and a discussion of our Initial concerns in 
the direction of inquiry. 

Underlying the RSE approach generally are the following: 

1. A valid state-variable model 

2. A choice of some specific criterion of optimality of estimation 

3. An objective of optimal estimation based on all observations 
from some Initial time through the present. 

The algorithms that result are characterized generally by both an evolving 
state estimate and an evolving covariance-of-estimatlon error. matrix 
(assuming model validity). The importance of the latter is underscored 
by the presence of the error covariance as a (matrix-) factor In the 
processing (Kalman) gain for new observations — as was the case In the 
loca My-optimum estimation algorithm (though the error covariance used 
there for gain calcul'ations was a static quantity), 

The updating of the error covariance in RSE, usually in a ’'down- 
ward*' direction as the data base grows, is a distinguishing feature of 
RSE and one that makes model validity so Important. For linear systems 
with a valid model, for example, the algorithm developed under the 
least mean squared error criterion produces the exact estimate value 
.(without approximation), and the estimator is finite dimensional and 
globally asymptotically convergent. If the mode! used is not sufficient- 
ly accurate, the expected and characteristic downward drift may be 
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present in The computed error covariance, but the latter may bear little 
relation to the covariance of the true error; a form of filter divergence 
results. 

When the system is nonlinear, the exact estimate in general is 
infinite-dimensional, necessitating for purposes of implementation some 
form of approximation, which, in turn may affect the modeling validity. 
Successful extensions of RSE to nonlinear systems In general are flnlte- 
dtmenstonal algorithmic approximations which are convergent In the 
locale of the true state value; from this standpoint RSE appears as a 
promising approach to tracking algorithm development for optimal MLS 
rece i vers . 

A valid model Is clearly a crucial factor in applications of RSE, 
and requires careful selection of many parameter values. To ease this 
task somewhat by reducing the number of parameters, one might be tempted 
in some cases to remove the stochastic forcing function from his model, 
i.e,, to model nis system as either unforced or forced by a knov^n 
deterministic function. The first implication of this is that the 
conditional expectation of the state at one point in time, given the 
state value at some other point in time, is not a random variable; the 
final consequences may be disastrous. With only observation noise 
in the model the calculated error covariance matrix tends toward singu- 
larity (actually may approach the zero matrix) as time and the (Implicit) 
data base grow arbitrarily large. The Kalman gain then becomes very 
small, limiting severely the Influence of later observation. In this 
situation the estimator is said to be in "data saturation" and is 
especially prone to the filter divergence syndrome, caused now, in 
part, by computational errors, such as round-off, etc. 

Two ways to avoid data saturation are as follows: 

1) Use static noise of sufficiant intensity in the model 

2) Use a limited memory approach to restrict the effective data 
base to the most recen"^ observations. 

The first approach Is the usuai one taken in RSE app 1 ications; the 
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application of locally optimum estimation to the MLS problem, as described 
previously. Is an example of the second approach. Both approaches 
effectively place a lower bound on the Kalman gain used In processing 
new observations. 


The digital computation of an evolving error covariance matrix 
presents certain problems also, arising from the f i nite-word length 
structure of digital machines. Covariance matrices in genera! are 
characterized by the symmetric and non-negative definite properties of 
matrices. Digital computational errors arising in a RSE algorithm can 
be categorized into two groups In relation to their effects on filter 
performance; 

1. Those that degrade the performance only slightly In terms of 
convergence and tracking properties. 

2. Those that Induce filter divergence. 


Computational errors made in the estimate extrapolation and update 
probably are of the first type and are remedied in subsequent updates. 
Computational errors made in the covariance matrix extrapolation and 
update tend to be of the second type; particularly those associated 
with the loss of the symmetric and non-negative definite properties 
Two ways to circumvent this problem are as follows: 

1. Integrate the covariance matrix equation off-line, if possible, 
on a high-precision machine and use the resulting steady-state 
solution as a matrix constant in the estimator. 

2. ' Reformulate the error covariance propagation problem In terms 

of propagation of a triangular ma+rix square root of the error 
covariance. 

The first approach may result in substantially suboptimal performance 
during certain periods of time, but provides additional benefits of 
computational simp I i city (on-line) and no possibility of data saturation. 

The second approach insures the (now implicit) covariance matrix will 
have the requisite matrix properties, but at a moderate i ncrease in 
computational complexity. Optimality of the estimate is maintained through- 
ouv all phases of its evolution, however, a benefit that must be weighed 
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carefutly considering the relative transience of multipath interference 
in the composite received signal. 

Finally, it is noted that stralght-forv/ard application of RSE 
theory to the MLS receiver problem gives algorithms which assume the 
carrier phase is known, as with the first application of the locally 
optimum estimation theory. Extension of the RSE theory is needed to 
treat the case of interest where the carrier phase is a parameter of 
no intrinsic interest, modeled as a uniformly distributed random 
variable on Extensions would be desirable also to the cases 

where the signal had been preprocessed with linear or quadratic detectors 
(recti fier- types), and consequently baseband observations were available 
but which were corrupted effectively by signal dependent noise that was 
neither Gaussian, additive nor zero mean, as discussed under locally 
optimum estimation. 

Conclusions 

Both the locally optimum estimation and recursive state estimation 
approaches possess both good and bad features in relation to the MLS 
receiver design problem. It Is expected that both algorithms should be 
considered for application in the final design. RSE could provide the 
desired extrapolation between scan periods and into fades, as required, 
given a valid state model. Locally optimum estimation might then be 
used to provide needed estimates of the state model parameters, which 
may themselves vary with time but knowledge is tacking of the associated 
laws of evolution and statistics. Much additional v/ork Is needed, of 
course, to obtain quantitative design and performance data on such 
a scheme and to describe the effects of such partitioning of the joint 
problem of identification-state estimation. 

Finally, the multipath propagation disturbance should be put back 
into the model and its effect on the estimation algorithm determined. 

The preferable type of algorithm would be an adaptive one. A non- 
adaptive approach vlould appear possible also In which the averaged 
performance of a class of multipath environments is optimized, but its 
performance in any one environment could be poor. The adaptive approach 
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is expected to be one in which the "gates" of the tracking receiver are 
modulated in some way by estimates of parameters associated with the 
multipath Interference model. 
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CHAPTER IV 

SUBOPT I MAL DESIGN AND SIMULATION EVALUATION 

A suboptima I receiver study was undertaken In order to consider 
some classical receiver designs involving envelope detection and 
subsequent processing by a form of early-late gating. This provided an 
opportunity to Improve the computer simulation as v/ell as to produce a 
candidate receiver design v/hich could be used as a performance standard 
in evaluating future receivers. No attempt was made to optimally process 
the envelope-detected signal. Simulation results, however, indicate 
that such a course may well be worthwhile. 

Receiver Algorithm Design 

Early design involved the development of a centroid receiver which 
utilized the relation between envelope pulse centroids and the A/C angular 
coordinate 0 . This receiver recomputed centroid positions (and thus 6) 
for each TO-FRO scan, rather than computing error in the estimate of 6 
from the previous scan; I.e., it was not a tracking receiver, and since 
It processed the signal envelope received over the entire scan interval, 
it had no multipath suppression ability. Several forms of early-iate 
gate tracking algorithms were then considered, with a square-gate version 
selected as the new receiver design. A development of the square-gate 
tracking receiver will now be presented. 

The square-gate receiver algorithm produces an estimate of the A/C 
coordinate by computing the error in the estimate from the previous 
scan. Let 0(t) = 9(t^) on the kth scan Interval (from 2-16), and let 
0(k) represent the receiver estimate of 8(t|^). 6 is formed as follows; 

0 (k + 1) = e(k) + A6(k) (4-1 ) 

A0(k) = /Jp[j(T)gCT,0(k)]dT (4-2) 

where t Is the local scan time at the A/C (from 2-23), Pq(t) is the 
direct-path signal envelope, on the (kt I )st scan, and gCT,0(k)H is the 
gating function (see Figure IV-I). The envelope functior is given by: 
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Pj^c-c) ^ apCe^(t) - 


(4-3) 


e,(x) = 


0Q + J2t 

0Q + fiTs 
6g + nT^ 


: 0 < T < T 

- s 

: T^< T < Tp 

s r 

Tp): Tp < T < T 


(4-4) 


where p<*) is the envelope selectivity function defined by (2-1), In 

Figure IV- I, tp| and tp^ are the envelope centroids. The receiver gates 

are centered at t,, and t_ which are estimates of tp, and tp- based on 

- y 1 92 I ‘ 2 

6(k). Each gate has height G and half-width w so that the gating function 


becomes : 


G; tg| - w < T < or tg 2 f T < tg^ + w 


gCtjOCk)] = \ -G; t < t < t + w, or t„ - w < t < t, 

1 Q I Q i Qo “ V 


0 otherw i se 


(4-5) 


t„ and t,^ are the values of t when 
g 1 9s A 


= e(k) and are given by: 


9Ck) - 0( 


(4-6) 




6ft - 6(k) 


(4-7) 


Define the error term 9 (k) as follows: 

e 


0 (k) = 0(t, ^,) - 0(k) 
e k+ 1 


(4-8) 


It follows from (4-1) that 0(k) ts an estimate of 6^(k). A relation 

e 

between the gate height and width Is desired which will drive C9(t. , ) - 0(k+l)I] 

A Kt I 

to zero by equating A0(k) and 9^(k). From (4-2) we have; 

0 


A9(k) = L * gCT,0(k)]p (x)dT + / , gCx,0(k)^p (T)dt 

Tg|-W u ^92“^ ^ 


(4-9) 
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(4-10) 


= s/!„Pd<^ + - ®/oPd''' • ®Od'’' + +92*''^ 

+ S/”pp(T + tg^ldx 

P|^(t + tgj) Is found by combining (4-6), (4-7), (4-4), and (4-3). Coupling 
this result with (4-8), we can rewrite (4-10) as: 


A6(k) = aGd/^ pCS)T - 6g)dt - j!!p(£3T - 0 )dT - p (-£2t - 0 )dt + /|!!fp(-J2T - 6 )dTl! 


(4-11) 

A 

Note that If 0_(k) = 0, A0(k) = 0, as expected. Assume 0_(k) Is small. 

© © 

Then: 


p02T - e„) = p(C!T) - 


(4-12) 


/°^p(inx - e^)dT : /°^pC±sit)dT + ^ = P(w) ? ^ [p(o) - ptEwl] 

(4-13) 

Substituting this result Into (4-11) yields: 
e_(k) 


A6Ck) = otGC (4p(J2w) - 4p(0))3 


-4aGe (k)Cl - p(flw)^ 
© 


(4-14) 

(4-15) 


Equating 8(k) and 6^(k), we have: 

© 


G = 




4aLl - p(S)w)J 


(4-16) 


It Is assumed that the envelope selectivity function p(0 will be known. 
Thus, given a knowledge of the signal amplitude a, the gate height G 

A 

can be computed for any desired width w such that 6(k) = 0„(k) (assuming 

© 

that 0 is sufficiently small). Although a is not known by the A/C, It 
can be estimated by noting in Figure IV-I that 


a=PpCtp^), i = i,2 


(4-17) 
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An estimate of a can be obtained by averaging signal envelope samples 
taken at tgj and tg^, which are estimates of tp^ and tp^. This estimate 
can then be substituted for a In (4-16) to compute G. 

The square gate receiver was developed without considering the 
effects of additive noise or multipath upon the received signal envelope. 
This approach can be defended In light of the following: 

1. The positive and negative areas of the gating function are 
equal. It Is therefore expected that the effect upon the estimate of 
noise occurring In a positive gate will be roughly cancelled by that of 
noise occurring In a negative gate. 

2. Any multipath distortions in the envelope occurring outside the 
gates will be Ignored. 

The receiver tested in the simulation can be described from (4-1) 
and (4"2) as: 


0 Ck +1) = e(k) + /gM(tj^^j ,t)gCT,0n<)3dT 


(4-18) 


where gC*) Is defined by (4-5) and Pq(t) has been replaced with the com- 
plete signal envelope MCt|^_i_l ,t), as presented In (2-53). On the kth scan 
a , is estimated by averaging ^ MCt^^tg^) (see 4-17), and G 

is then updated; 


G ^ 




2LM(t. ,t ) t M(t. ,t )J X Ll - pCRW)J 
K g I ‘^92 


(4-19) 


It was decided to set w equal to the time between centroid and fi.st zero 
of an envelope pulse of Pq<t). This selection of gate width seemed in- 
tuitively good In that the gate would include most of the direct signal 
envelope and yet be narrow enough to exclude most multipath. 


The computer simulation involved a discrete- time version of the 
local scan time t, so that the square-gate receiver operated on samples 
of the envelope M(tj^,x). This operation will be described In the next 
section of this chapter. 
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simulation Model 1 ng 


The simulation objective was to evaluate envelope-detector receivers, 
making it necessary to mode) the l.F. signal envelope M(t^,T) In (2-53). 
The signal modeling was baslcly that of Chapter II, with the following 
changes: 

I. The ^'received scan interval*' of (2-22) became 
r(t, ) 

t, + — < t < t, + ^ + T (4-20) 

k c - - l< c 

r ( t ) 

in which fairly insignificant — terms were dropped, in terms 
of local scan time x this interval became (from 2-24) 

0 < T < T (4-20A) 


tp^, defined in (2-25), likewise became 
r{t, } 

to = t, + ^ (4-21 ) 

«K k c 


2 . 


3. 

4. 


The signal amplitudes a and a. in (2-42) and (2-44) were 
assumed constant. This is valid If the A/C and reflector 
ranges do not vary appreciably. 

tj^ was dropped from the argument of in (2-42) and (2-44), 

as 0^ is a periodic function which reinitiates at t = tj^. 

The range and angle coordinates of the A/C and reflection points 

• « « • 

were assumed to be linear functions of time; i.e,, r, 0, fpjj 


and were made constants. 


(4-22) 


The local scan time t was discretized to simulate sampling of the 
received signal envelope M(t|^,T). Starting at x = 0, the envelope was 
sampled across the entire scan interval at a rate equal to the IF band- 
width B|p(160KHZ). From (4-20A) the number of samples taken in one 
scan interval would be 


J j - T X B| p + I 

The local scan time at the jth sample would thus be: 


(4-23) 


4 ! 


T(j) = Cj - 1)6 


(4-24) 




where 6 = l/Bjp. The envelope function can now be expressed in 

terms of the discrete local scan time by substituting (4-23) Into 
equations (2-47) through (2-50) and incorporating the new changes In 
modeling enumerated above; 


0 ) r 
c 


Yp (tj^,j) = apCe^((J - !)6) - 0Ct(^):cosCe(t(^) - (j - 1)6] (4-25) 

A 


(4-26) 


= Pp(t|^J)cosH'p(t(^,J) 

Yq = “pC0y^(CJ ' 1)5) - e(t^)]sin[:B(tj^) - (j - |)6] (4-27) 


“c*- 


- PDCtj^»j)sinWp(tj^,j) 


(4-28) 

Ar. (t ) tu r^(f ) 

Vr (t,,J) = I a.pC.fU - 1)6 - _!^) - e (t,)3cosi:e,(+^) - 

^ ' (4-29) 


I Pr (t^,J)cos'l'R^(t^^,j) 


i k' 

Ari(tk) 


yrs<vj^ = I«ip[:0Aca - 1)6 — ^^) - 


- IPRt 


(4-30) 


w rj.(t.) 

eR,(t^)]stnD,(t^) - (J 

^ (4-31) 


(4-32) 


Note that the RF phase terms and H'r. Ct^pj) are linear functions 

of time and can be expressed: 




B(tk>: j = I 


(4-33) 


'i'oCtk^j-D + A'Pp-. 2' < j < 


Where ATp, = - — 6 
□ c 


(4-34) 


j - 1 ) 5 ] 


-1 ) 6 ] 
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v<VJ> 


r j = ' 

L'*'r ‘VJ’” "R|‘V= 2 < j < 

“VTi<+k> 

where A'I’d (t,,) = — — 6 

Kj k c 


(4-35) 


(4-36) 


The envelope M(t^,j) now follows directly from (2-54) through (2-56); 


c c s s 


(4-37) 


where n (j) and n (j) are the quadrature noise components of (2-51). 

*-F< W 

The simulation must compute M(tj^,j), j = 1 to Jj, for every TO-FRO 
scan of Interest. This also involves computing A/C and reflector coordi- 
nates and velocities at t^ (and in some cases t^j^), as these appear as 
parameters In the equations Just developed. A description of the 
simulation for the kth scan will now be presented. 

Let t, (k = I , 2, . . ,k ) be evenly spaced, and define: 

K msx 


1*1 in = *^1 ■ ■*'i, I "75 milliseconds) 
UD K k- 1 


From (4-22) we can update the A/C and reflector coordinates: 

rCV = "^Vl’ ^ "‘■^UD 
0CV -L0*T^, 

4>(\) = i-Tuo 

^R. ^'*'k^ '‘Rj‘ ^UD 

®R, ®R,‘ ^UD 

" ‘^’R^'^k-l^ ^R,* 

I i I 


(4-39) 
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For each reflector we compute r^ (t ) from (2-12) through (2-14): 

I I K 


= coS(}iCtj^)coS(i>Rj (t|^)cosCe(t|^) - sin(|>(tj^)slniijp. Ctj^) 

(4-40) 


'•AR,«k> = Cr^'V - 'i‘+k>'-‘V^R/V ^ 
•■t '+k> ' '■r, '■ar, '+k> 


4r,(t|^) = - rCt,^) 


(4-4 1 ) 
(4-42) 
(4-43) 


The time derivative of (4-41) can be shown to be; 


r(t, ) 

r (t ) = ^ 

ARj K ‘‘AR^’^’k 


•"r, 


) ■*■ r^„ (t,J ^""R, " 


ARj 'k' i 


(4-44) 


"t, ^-^k^ = "R. "AR, 


(4-45) 


From assumptions (2-20) and (4-21) we obtain; 


r(t ) 

r(tR^) = r<t,) + ^ 

r(t, ) 

'"t, ‘■^rk’ ' ■‘t, <+k> -T- 


^"i ^"**RK^ " "t, ^’’’rK^ " 


(4-46) 

(4-47) 

(4-48) 


(t,^) Is computed from (4-36)., Set H’,^(t,^, I ) to 3(t,^), a uniform 
random variable on C-ir^uH supplied by a random number generator. From 
(2-40) and (4-35) compute: 




*R,'V> = = ®'V 

(4-49) 


44 


Now can be computed for each sample on scan K. 

For j = f to Jji 


= ap[0^{(j-i )5) - et+kO 

' 

C4-50) 



(4-51) 



(4-52) 



(4-53) 

Ar (t ) 

= a,pCe^((J-l)S ) 

- ®r/V^ 

(4-54) 

C I 1 ‘ 


(4-55) 

\ ° iPR/VJ’^'^Ri'VJ’ 


(4-56) 

■’’Ri'VJ' ^«R,<V 


(4-57) 


Obtain quadrature noise samples n (J) and n Cj) from Gaussian random 

w 5 

number generator. Compute M(tj^,j) using C4-37). 

In the FORTRAN computer simulation, the antenna scanning function 
e^(*) Is computed by the function subprogram THA(T), where the argument 
T is computed in the main program, 9^^*) defined by C4-4) and illus- 
trated in Figure IV- I. The angular selectivity function p(0) is com- 
puted by the function subprogram P(THETA), where THETA is computed in the 
main program. is defined as follows, from reference D, pp. 2-213, 

2-214]: 


pC0) 


where a = 2,46 C4-58) 


^ cos( ^ a) 






4 


3 I 


a - 


V- 
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The simulation program computes M(t|^,J) for J=l to Jj and stores 
these samples in array M. M Is then passed to subroutine RCVR, which 
contains the square gate tracking receiver algorithm. The algorithm 
updates the estimate by computing the error In the estimate of the 
previous scan. It differs from the receiver developed earlier In this 
chapter as fol lows: 

1. The signal envelope now exists in discrete form so that the 
integration of (4-18) is performed using a trapezoidal approx- 
imation rule. 

2. Discontinuities In the gating function (see Figure IV-I) have 
been forced to occur at the sample times, so that the gates 

are well defined. Otherwise, the receiver would produce the same 

A 

value of A6(k) In C4-I8) regardless of where between two samples 
the gate center occurred. 

RCVR operates on the samples of array M one by one, with the result 
stored in array EST. EST thus contains an evolving estimate, where ESTCJy) 

A i 

represents the final estimate 0(k) for the kth scan. 

After computing 6(k), the simulation program advances to the next 
scan and recomputes the envelope and estimate. This is continued for 
as many scans as desired. For efficiency in computing the simulation 
has been carried out in two FORTRAN programs. The first, MLSRCVR, 
performs all simulation calculations and stores results on file. The 
second program, MLSPLOT, plots these results using a GDC COMPLOT DP-7 
digital plotter. Emphasis is placed here on MLSRCVR. The following 
plots are made: 

1. A long-term plot, which plots for every scan K In the simulation: 

A. the composite signal envelope at the time of direct-signal 
centroids 

fs 

B. the error in estimate 6Ck) 

2. A plot of signal envelope M(t|^,j) as a function of local scan 
time for six selected scans. 
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3. A short-term plot, which for one selected scan plots the 
received envelope M(tj^,J) as well as the evolving estimate 
as functions of local scan time. Also plotted is the evolving 
estimate which would result from a received signal containing 
no noise or multipath, 

A flow chart for program MLSRCVR Is presented In Figures IV-2 through 
IV-3, while a basic flowchart for MLSPLOT appears In Figure IV-4, Pro- 
gram listings are contained in Appendix A. It should be noted that 
variable names used In the flowcharts and programs often differ from 
those used hare. 


Performance Evaluation 

The square-gate tracking receiver was tested under the following 
conditions; The A/C with coordinates r = 10 K.ml., 0 = 30®, 4» = 3® 
was approaching the runway with airspeed 300 knots (0 = | = 0). A 
single specular reflector with coordinates rp^^ = I.O N.mi., 6p| = 53®, 
cj)pj| = 1.85®, v/as following the circular path rf^| = 0, 0pj| = -3®/second, 

=5 0. One second later (scan !4), the reflector was at 0pj| = 50°, 
so that the signal envelope at the A/C due to the reflector was Inside 
the receiver gates and coincident with the envelope due to the direct 
signal. After another second (k=27), the reflector was at 6p| = 27° 
and the signal envelope due to reflector was outside the receiver gate 
again. This was considered to be a realistic test for the square-gate 
receiver In that in-gate multipath interference intuitively represented 
a worse case. The simulation was run for 27 scans with the above initial 
conditions, with signal amplitude ot = 1.0, and multipath amplitude ~ 0.8 
and for several values of Initial RF phase difference C3j(l) - 8(1)1]. 

Also made were two trials with no reflectors and differing signal-to- 
noise ratios. 


All simulation trials were made on the CDC 6400 computer. Plots 
of receiver performance may be found in Appendix B. Results are tabu- 
lated in Figure IV-5. Note that with a 20 db signal -to-no I se ratio 
and no multipath the receiver was able to estimate 6 with an R.M.S. error 
of 0.015®. R.M.S. errors are shown for the trials with multipath, yet 
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Input Data & Set Constants 


Input Initial Conditions for A/C S, Reflectors 
Input Parameters for Short-Term Plotting Decisions 



Update Signal to Start of kth Scan 

Compute Baseband Signal As It Evolves over kth 
Scan Interval and Store Samples 


Compute Receiver Response As It Evolves over the kth Scan 
Interval and Store Final Results in Long-Term Plotting Arrays 


if Desired, Store Short-Term Plotting Arrays on File 


Advance Geometry to Start of Next Scan Interval 


Compute RN!5 Error in Angle Estimate for k Scans 

^ max 

Store Long-Term Plotting Arrays on Fi te 




Figure lV-2 Flowchart for Simulation Program MLSRCVR 

(Note: Details of blocks identified with encircled 

letters are given on subsequent pages.) 
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Figure iV-2c Updating Signal to kth Scan 
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Enter Subroutine RCVR (Receiver Algorithm) 


Input: 0^g^(Flnai Estimate from Previous Scan) 

M (Array of Signal Envelope Samples) 


Output; EST( Evolving Estimate Array) 

(EST(Jy) = Final Estimate of 9) 


Fi I 1 Long-Term Plotting Arrays 


TMK(k) = 

CS(K) - MlN(CSi,CS2) 


THER(K) = 6 - EST(J^) 


^RMS “ ErMS CTHERfK)]' 


Figure IV-2e Estimate Calculation 


53 
















55 






Advance A/C Geometry to Start of Next 
Scaa Interval 


■■tk ' ^tk ''^UD 


e =0 + 91 . 


(j) = <j( + ^T, 


Advance Geometry of Ith Specular Reflector 
to Start of Next Scan Interval 




6r(i) = 9^(1) + e^n)T^^ 


♦r'I’ ”■ ♦r*" ♦r‘”'^ud 



Exit From 
Block G 


Figure IV-2g Advancing Geometry to Next Scan 


55 







Input Array tl 


OEIGINAL PAGE IS 
OP POOR QUALET^ 


Romark.; lj_, T(, Oo, Hi 

t, fi. Stnrad IntornBUy 

t 

AO = tt = 0 

Jp = tFlXC(6 - OoWtilT) t 1,5] 

a = mtjp - 1) + 8o 
S| = f‘ 


1.-1 rJ j 

Jo = Jp - '0 
Jl “ Jp + 10 


J. - Jo + 1 

-H 


40 s 40 + l/Zig,xCM(j- I) + M(J)D 

■ esTij> * i ♦ 40 . 



Jp .= IFIXCCTp + (6j ' + 1.53 



g = pD - <z.4ae)®3/(ZDi 


6 = -n/(4g) ■ ■ 

' t ^ 

r sw j 

Figure lV-3 Square-Ga+e Receiver Flowchart 




























Figure 

Jobname 

SNRIdb.) 

No, 

Ref 1 ectors 

RF Phase 
Difference at 
Coincidence 

R.M.S, 

Error 

Comments 

B-l 

MLSRCMP 

20 

0 

- 

0.015® 


B-2 

MLSRCAZ 

8 

0 

- 

0.059° 


B-3 

MLSRCJO 

8 

! 

180® 

0.247® 


B-4 

MLSRCFG 

20 

1 

o 

O 

00 

0.203° 

Fade; Greatest Error 
Coining out of Fade 

B-5 

MLSRC5C 

20 

I 

0° 

0.162° 

Phase Enhancement 

B-6 

MLSRC8E 

20 

I 

0 

o 

0,173° 


B-7 

MLSRC70 

20 

1 

270® 

0.201° 



Figure lV-5 Simulation Results 
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those quantities may not be very useful comparison measures in genera!. 
This Is especially true for differing multipath environments, since 
R.M.S. is a time-averaging process and the muitipath Interference may be 
non-stat ionary, A more useful measure may be peak error, which is 
avai I able on the plots. 

Note in Figure B-4.a how the signal fades when the direct and 
reflected signals are beam coincident but 180° out of phase (scan 14). 
Also note that the greatest error in the estimate 0Ck) occurs not at 
the fade, but as the signal comes out of the fade, in Figure B-4,b 
one can see how the reflected signal envelope moves across the direct 
signal envelope as the simulation progresses, and how at scan 13 (I scan 
before coincidence) the two signals have nearly cancelled each other. 
This is contrasted with Figure B-5.b, where phase enhancement causes 
the signal envelope to increase In amplitude on scan 13. Also note 
that in general the simulation trials with multipath have very low 
error near the beginning and end of the trials. This is expected, since 
a reflector removed 3° in azimuth from the A/C would produce an out- 
of-gate signal envelope pulse. 

Conclusion 

This study has produced and demonstrated a receiver algorithm of 
basically intuitive design which has some commendable performance 
features. The value of the simulation developed for evaluation of 
candidate algorithms was also demonstrated. 

The study of this algorithm suggested two options for further 
consideration; 

1. An algorithm which optimal ly uses the detected envelope to 
estimate the A/C angular coordinate. 

2. An algorithm not closely coupled to the exact form of the 
selectivity function p(0, and hence. Inherently suboptima I, 
but perhaps less sensitive to variations In the shape of 
p(0 arising from cross-coupling of the angle channels and 
slte-to-site design, and installation variations, etc. 
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The simulation appears to be validated In the results presented. 

The principal use in the near future of this tool will be In the study 
of optimal tracking algorithms. Some refinements will have to be made 
as we consider both quadrature detection and envelope detection 
estimation algorithms. Also, a linear evolution of A/C and/or reflector 
coordinates seems more plausible in a Cartesian frame than in a 
spericai system. This modification will be made to test the validity 
of assumption, equations (2-15) through (2-20), made in signal modeling. 
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CHAPTER V 


CONCLUSIONS AND RECOMMENDATIONS 

A summary has been presented of both theoretical and simulation 
study results, including 

1. The development of detailed signal models cast In the time 
frame of the A/C receiver; 

2. Preiiminary optimal estimation studies. Involving first an 
application of the locally optimum estimating algorithm to the 
direct-signal and noise specialization, then a critical over- 
view of some relevant characteristics of recursive state 
estimation approaches; 

3. Design of a suboptimal square-gate tracking algorithm which 
accepts the baseband output signal of a standard linear AM 
demodulator, and then evaluation of this algorithm in simulation. 

The greater significance of the work done to dale is in identifying 
and bringing Into focus concrete problems end problem areas that should 
be -onsidered, and In suggesting methods "of attack. 

First among these is the modeling of the total reflection component 
of the signal; the modeling results described gives support to the 
belief a random field defined with respect to a low-dimension parameter 
space may be a suitable model. This would motivate a multipath-adaptive 
approach. 

Another problem stems from our ignorance of the law of evolution of 
the A/C angular coordinate and the general and unknown variability 
with time of the law. One approach which should be considered involves 
simultaneous model Identification and state estimation, possibly in a 
"layered" algorithm with different criteria of estimation optimality 
for parameters on different levels. 

The study of processors using the output of AM detection should be 
broadened to Include both linear and quadratic AM detection, and also 
an objective to process optimal ly the detector outputs in calculating 
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an estimate of the A/C angular coordinate. The most attractive estimation 
algorithm for the MLS receiver may well come from such an approach. 

Final ly, in support of the prototype construction effort to be done 
In the project continuation a number of decisions must be considered 
and made, dealing with the choice of algorithm, choice of microprocessor 
and choice of interface equipment. Synchronization and acquisition 
details must be clarified and a system plan worked out which wl I ! allow 
the microprocessor to perform all its jobs in the requisite real time 
frame. These tasks are under consideration presently. 
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APPENDIX A 


COMPUTER PROGRAMS 
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I 


1 


1 


J 



) 


PROGRApi MLSRCVR(IMPUTiOUTPUT,TAPE20) 

C 

c 


■ 

c 

c 

MICROWAVE landing SYSTEM RECEIVER S.H, IRWINiJR. 

REVISED 10/3/75 

9/1/75 

c 

c 

MLS RECEIVER SIMULATION PROGRAM, IN THIS PROGRAM THE RECEIVED 

signal Is 

c 

c 

generated and recCivPTperformance evaluated, output data is 

PERMANENT FILE, WHICH Is THEN ACCESSED BY A PLOTTING PROGRAM. 

stored ON A 

^ __ _ ^ __ __ . 1 

c 

c 

THE following MUsT BE ON DATA CARDS AJ THE END OF tHE PROGRAM! 


c 

c 

lOPT 11 

NAME(I)*I = H3 3A10 


c 

c 

iNRcVRCI} tl = l*3) , (NRCVR2( I > t I=1»3J 6A10 

CHAW 3A10 

. 

c 

c 

ERSPEC E20.1 

FC,TUDtTHO,lHS,TS,TR,BEAHWTH 7E10.2 


c 

c 

BIFtALPHA 2F10.3 

ALPHAM FID. 3 


c 

c 

SNR F10.3 

THES E20.1 

-■ 

• 

' 

c 

c 

(RAtl) ,1=1,5) 5E15.1 

IMAXtKMAX 2110 


c 

c 

RTK,V,THETA»OTHtPHI,DPH 6(E10»1) 

RRTK ( D , VR ( 1 ) , THETAR ( I» , OTHR ( I ) , PHIR ( I ) , DPHR 1 1 ) , ALPHAR ( I ) , DI FO { I) 6^1 0 . 1 

c 

c 

(FOR 1=1 TO IMAX> 
{PL(D ,1=1,5) 5E15.1 


c 

c 

c 

LOGICAL PlTSCAN,FLAG 


Real MtnPi\iHtHc»fisiNc*ws 

DIMENSIOW M(2000)»OM(2000)*RM(2000)tTM(2000 ) lESTCEOOO ) »E:sI(2000 » 

DIMENSION RRTK(10»tTHETAR(lo)7ErTHR(lQ)iPHlR{lo)iDPHR(lO) 
dimension AlPHAR(ID), OIFO(IO) ,TDR(10) ,PHASER(10) ,DPHASER(101 

DIMENSION VR(10)»NRCVR(3) iNAME( 5J fFlTO(iO) 
dimension CS(30)fTHK(30JfTHER(30} tTA(2000) 

dimension RA( 5) *RA1 (5> ♦RA2(51 »NRCVR2(3) ,PL(5) 
COMMON TRiTSiTF»THotTHs*OMGiTAU 

COMMON/PLf/KMAX 1 ERR t PL 
common/sel/beamwth 

c 

COMMON/TRACK/THES, RA 


c 

c 

INPUT Plotting option iopt, where; 


c 

c 

I0PT=1! generate ALL »LOTS 

I0PT=2! GENERATE LONG-TERM PLOT ONLY 

. ^ 

, 

c 

c 

I0PT=3: GENERATE SHORT-TERM PLOT ONLY 

* 

. tr; 

READ ll,IOPT 

11 FORMAT (ID i 

C ' ' i 

C 

c 

c 

INPUT TITLE INFORMATION FOR PLOTTING PROGRAM 

1 
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PBMJEDING PAGjS BLAJ^E NOT FILMEID 


ORIGINAL' -- 
OF POOR QU A-IiTi 







n n o n n o n o n n n o o |o op n r* n jo op n 1 npo,nnpn 



NAME=UsCRS NA ME 

NRCVRaWAME OF RECEIVER ALGORITHM 


' j R C y R 2 = A n D I T I p N A L . R C V K_ INFO R H A T 1 0 W ( J f_ P > 

CHAN=MLS channel OF INTEREST 
EJLSEEC=EM£lfl_ae-EC 

READ lOt (NAHE(J> 1 J=l»5> 

in FORHAT(bAlO) 

REJp_ 10 _i t N RCVR LiiL* J=l«3) t <NRCVR2 (n *1=1 1 5 ) 

READ lb t CHAN 

R E AP_ 2 0 j, E R SPEC 

2*ri FORM AT (E2 0,1) 

JB=JObiNAME(X) 

DAY=DAtE(X> 

_ WRIT EJ^ Pi ^ pPTj N A^H Ej, N RH/ R»NR C V^2 t CHANtE RS PECi JBi DAY 


» **************** *******=t^ *********** 
* . * 
* * 
■* initialization- input DATA, SET CONSTANTS * 
± * 



!»:;<;» **;»*****;^ ************************** *********** 


SET CONSTANTS AMD CONVERSION FACTORS 



csspeed of light in meters per second 

?1P;\tM=METERS PER NAUTICAL MILE 
S PH=SECONDS PER HOUR 

C=3.QE6 

PI=3,1A15926536 

MPMM=iaS2. 

SPH=3600, 


INPUT PROPERTIES OF MLS ANGLE FUNCTION 



FC=CARRIER FREQUENCY (HZ) 

TUOtUPDATE TIME) -TIME BETWEEN rKE^TAI 

the next scan of the same angle function. 

THO ( THETA-0) -SCAN ANGLE AT TNE BEGINNING OF THE TO SCAN 

THS(TH£TA-S)-SCAN angle at the end of the TO SCAN 

TSCSCAN TIME) -length IN SECONDS OF THE TO sCaN 

trcreference time)-time in seconds between zero-degree crossings 
OF the scan angle 

deamwth=beamwidth OF the scanning antenna pattern 

read InOiFCtTUl'j 'HotTHSi TSiTRiBEAMWTH 


IQO FDRMAT(7E10,2) 



- 1 




\ 


• 



c 

t 

initialize PLTsCAN and GAUSS 



c 

PLTSCAN=.F. 



c 

ZZ=GAUSS(SQRT(g, ) ) 



c 

c 

c 

c 

INPUT RECEIVER CONSTANTS 



c 

c 

BIF=IF BANDWIDTH 



c 

L 

ALPHA=PE'AK DIRECt SIGNAL STRENGTH 
ALPhAH=AMPLITUDE OF scattered MULTIPATH 



C 

C 

SNR=SIGNAL to noise RATIO IN DB 



READ RoO'lJXt- 1 ALPHA 
200 FORMAT(2F10.3) 

READ 300*ALPHAM 
READ 300tSNR 

300 FORMAT(F10,3) 

SIGMA=ALPHA/SQRT(2.o)*10.0**(-SNR/20*0> 

c 

c 

C 

c 

CALCULATfc AU60RITHH CONSTANTS PROM INPUT DATA 



c 

c 

tau=timf. in seconds between samples 


■ 

c 

c 

OMG( OMEGA ) =SC aN RAYE ( DEGREESZsE^ 

tf:;Time in seconds from start of to scan td rtart of 

Fro scan 

■ 

c 

c 

jt=total number of samples 



TAU=l./BiF 
OMG-( THS-THO )/TS 

C 

TF=TS+TR"2.*THS/0MG 

ti=sscan interval from start of to Scam to enq of fro 

SCAN 


c 

ti=ts+tf 



JT=l + IFXX(Tl!t^BIF+0,5) 
WCC=2»*PI+FC/C 

c 

C 

C 

c 

INPUT tracking RCVR CONSTANTS: 




c 

c 


c THES=ESTIMAfE OF THE A/C ANliULAR' POSITION C A L CU L ATED~ 0 M t HrTffCV I 

C RA= ARRAY' OF CONSTANTS PECULIAR TO A SPECIFIC TRACKING ALGORITHM (COMPUTED 
C 'EXTERNALLY), RA MUST BE COMPLE'TELY FILLED. IF 5 PARAMEfCKS ARE NUf 

C REQUIRED. LEFT-OVER POSITIONS ARE sET TO ?ERO. 

C 

C 

READ 53GiT)iES 

550 FORMAT(2E20.1) 

READ 340 * ( R ATiTTI s 1 . 5 ) 

340 FORHaT(5E15.1 J 
DO 350 I=l»5 
RAK DsRAt n 
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350 CONTINUE 

£ 

C 

C INPUT N UMBE R OF SCAN S AND RE FLEICTORs 

C 

c''Ki^AXrNLJ(iBER OF SCAN^ 

C IMAXsNUHDEP. OF REFLECTORS 

C 

READ 4D0< IMAXtK H AX 

400 FORMAT (21 10) 

WRITE(201 KMAX< IMAX i A lPHAi AlPHAMi SNR t SIBMA » Tl t J T 

C 

C 

C 

c input aircraft initial, conditioivs 

C 

_c ^ 

C RTK=RAWGE in nautical miles 

: V=AIRSPECD IN kno ts 

'c theta=a/c azimuth angle 

C DTHrT iME CHANGE IN THETA 

“C PHI=A/C elevation ANGLE 

C DPH= TIME CHANGE IN PHI 

C 

READ 500 tRTKiVt THETA fPTHt PHI iDPH 

500 format (6E10.1) 

WRITE(20) RTK. ViTHETAtOTHtPHI tDPH 

PRINT iSOlf {NRCVR( I) 1 1 = 1*3*) 1 (NRCVR2( I) t I = li3) 

1501 F O RMATt iHi* i 0X*3Al Dt5X*5Al0 ) 

PRINT 150 2tSNR*"fHETA.THES 

1502 FQRMAT(lH0,10Xt*SNR=*«F7.2<5X,*THETA=**F6«3t5X,it=QRlGlNAL ESTIMATE^ 
1+ 1 FO * 3 ) 

PRINT i503*KMaX* IMaX 

3.503 FORMAT (lhOtlOX**Nb, SCANS=*T1 t *5X * *NO, REFLECTORS=* • 2*2 ) 

PRINT l5Q4iI0PTiCHAN t ERSPL C 

1504 FORMAT«lH0,//i llXi+IOPT=*, I2*5x*A10*+CHANNEL*,5Xi*LRSpEC=+tElI,3) 

PRINT 1505* (NAME( I) *1 = 1*5) *JB. DAY 

1505 FORHAT(1HO*1DX**RUN bY * * 3A10 * 5X, ^-JOyNAViE . :c, Alo • 5Xi *dATES **A10) 

PRINT 1506»ALPHA, ALPHAM*BIF 

1506 FORMAT(1HO*10X**ALPHA= * i F6t 2 1 5X < *ALPHAM= F6. 2* 5X * *I tJArH'^YTlTf 

1H= *iF7.0** HZ*) ^ 

PRINT 1507 

15Q7 format ( lH0t//i20X**MLS FUNCTION PROPERTIES;*) 

PRINT 1508*FC*TUD*THo*THS*TR*TS*BEAMWTH 
15 OS FORMAT ( IHO t lOXi *FC=* , Ell * 3 1 3X , *TUD==(s i F6. 3 , 3X , *THO = l: ,F6 ,1 , 3X , +THS=# 
ltFS;i*3X**TR=»i=iEll,3*5X**rs=**EU,3t3X*lcBEAMWlbTH=*fF5t2) 

PRINT 1509 

1509 FOR'MAT (lHo*//*20X* ^iTRACKINGHTECEIVER CONSTANTS;*) 

PRINT 15lO«tRA(I)*I = l*5) 

1510 FORMAT(lHOtlOX**RA= ( *« 5 (2X *£11 . 3 ) * * ) * * ) 

PRINT 1511 _ 

1511 FORMAT (iH0"t7/72“0X»*A/C INITIAL CONDITIONS:*) 

PRINT 1512*RTK* VtTHETA*DTH*PHI*DPH 

1512 forma fTlHo * ldX"**R=** F6 , 2 * S'X* *V = *tF7 »2 * 5X f ^THeTT As* *F7 , 2 * 5X « *uTFr=*TF' 

16.2<5Xi*PHl=*iF6.2f 5Xt*DPHs*tP6«2l 

RTK=RTK*MPHM 

V=V*MPNM/SPH 
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1 


u a u'u u u o u <J<J a C3^ u u o u 


t)PHASE:o=:WCC>!‘V*TAU 


INPUT initial conditions OF THE'TTH MFLECTOR 


I F< IKAX.EQ«0) go to 1030 

DO 1000 IsldlHAX 

READ 600 1 RRTK ( I ) ♦ VR 1 1 ) ,THETAR (I ) , DTHR ( I ) i PHIR ( I » i DP hR ( I ) i ALPrAR { I ) 

rroiFom 

e OQ FORMAT( aClO.l) 

PKINT i513tl 

1513 FGRMAT(lH0*//i20Xf*INITIAL CONDITIONS FOR REFLECTOR* . 12 ) 

PRINT 1512 1 RRTK ( 1 ) , VR ( I ) » THET AR ( I) , DTHR ( 1 ) i PHIR ( IJ t DPHR { I > 

PRINT 1514,ALPHAR( I ) ,DIF0(I1 

1514 FOR7iAT(1HO,1C!Xi*ALPHA=*,F 7;2,10X»*PHASE DiF. FROM DIRECT SISNALtRA 

10. )_=^iE_l6.8) 

IFU.LE.3) WKITE(20) RRTK ( i ) , VR { I M THE TaR ( I ) * DTt:]R (I ) * PHIR { I ) i DPHRi 

11^ . ALPHAR ( I ) 

RRTKa 1=RRTK(I )*MPMM 

VH { I )J=VR ( I)*MPN M/SPH 

1000 continue' 

I F I IHAJ( . 3_)_ GO TO 1050 

i030' IJ=3-IMAX 

DO 1040 1=1, IJ 

WRITEtaO) C*CiCiC<CtC*C 

1040 cont inue 


input sHORt-tERM plotting constants? 


ARRAR' PL contains' ConST ANTS“ToR ACCESS BT“sUbR^ T I nE^ SHORT * W'hlcH'DEtrRMlivl'ES 

what scan will HaME a detailed short-term plot of receiver performance. 

IF ARRAY PL IS NOT NEEDED, if MUST 8E FILLED WUH ZEROs, 


1050 READ 340. (PLC I > . I=1.5> 

PRINT 1515. (PLTI) .1=1.51 
1515 FORMAT{1HO»10X»*PL= ( *. 5 { 2X. E16 ,6 ) . * 1* t ) 
lNC=KMAX/6 

1NC2=1 

ERHS=0,0 





TM(J)=(J-1>*TAU 
WxfTEtaO) TM 

DO bOOOO K=1(KMAX _ _ _ 

TKsCK-l)*TUD . . „ 

IU.U,=RTK^.C 

riK=TK+TDD 

Ft-RrK-»V*T00 

PHASED=:2.D'»'Pl*{RANr<0.0>-.5) 

PRIAi T 15l7' K.RtTHET A tPHI 

1517 FOKkAT(lH0t///*5X»*SCAN**r3*l0Xi>t'R(METERS)==t‘tEl6“,B,5X*»THE:TArijEG) 
l=*»Ei&,8*5X»*pHI (DEGJ =*tE16.B) 
lF(IMAX,EQ,d) *60 TO 2050 

00 2000 1-lflMAX 

ZETA=CGS<PHI)*COs(PHIR(I) )*COS<THETA-THe:TAR(I) ) +SINCPHI )*SIN(PHIR( 

ID) 

RAR=SQRT ( RTK»*2-2 . *ZETA + RTk*HRTK ( I ) + ( RKTK (I M **2 ) ™ 

VAR={RTK* (V-Z£TA*VR<I ) )+RHTKtI)*(VRtI)-ZETA»»V) >/RAR 
TDR{I)=(RRTK< I j+RAR)/C 

RR=RRTK(IUVFUI)*TDR{I) 

RAR=HAR + VAR+.TDHCI ) 

RT=KK+MAK 

DPHASER(U=WCC*IVR{ n+VAR)*TAU 
IFtK.WE.lJ 60 TO 1900 

DlF“DIFO(n ■ “ 

RO=R 

RTO(I)=RT ™’ ~~ 

GO TO 1950 

1900 DIPsOIFOtn-WCC^CRT-RTOdJ-H+RO) 

00 1800 OK=l,lO 

IF(AB*S(D‘IF) ,LT.2.0*PI) GO TO 1950 
TF(DIF.GT.O.O) DIF=DIF-2, 0*PI 

IF(OIF.LT.O.O) DIF=DIF+2»0*PI ~~ ^ 


laoo CONT INUE 


1950 

PHASER(I)spHAS£D+biF 

OlFO=DIFflS0,0/PI 


1518 

print 1518,I,RRiTHETAR{I) iPHIRUJ 

F0RMAT{lH0,10Xi*REFLECTOR+,I3,5X,+Rs:*i£lG.6,5X,*THETA = i4:,El6.6*5X«* 


iPHi='»i.Ei&,a) 
PRIWT 1519, DIFD 

1519 

FORMAT ( 1MO,10X,+RF PHASE DIFFERENCE FROM DIRECT SIGNAL-*, Eld, ti,*DE 
IGREES*) 


2000 
C 

CONTINUE 


C 

C GET sample numbers CORRESPONDING TO DIRECT PATH CEmTROIDS 

C 

2050 

JP1=1+IFIX( (THETA-THO)/{TAU*OMG) +D,5) 


C 

JP2-1+IFIX( { (THS-tHETA)/bMGTTF)/TAUt6,S) 


C 

C 

c 

c . 

* « 


c 

c 

* * 
St; COMPUTE baseband SIGNAL FOR THE KTH SCAN AND STORE SAMPLES * 


c 

c 
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i 

V 


C ♦»i» + **'f^’(:’t:t>t;#!(:*******!<i*»’(f + **’t<***** + + ***************:(:**-*;|i*** ******* 

C ■ 

c 

DO 10000 J=lfJT 

' Tsi=rj-niTAU ‘ ‘ ‘ ‘ ^ ** ■ 

T = T1K + TSI _ _ ' 

TA('J)=THA(Tsi ) 

F^D =P(TA < Jl-THETA) 

YDC = PD+cbs(PNAS£:D) 

YDS=pU*SlN{PHaSEO) 

■'THASEDiPHAS'EWDPHASEb ■ 

ISC£.0. 

YSS=0.’ 

YMC=YSC 

YMSsYSS 

IFUMAX.EQ.O) go to 3100 
DO 30'0 0" I=iilf-1AX 

PR=ALPHAP( n*p(THA tTsI+TDD-TOR(I > )-THETAR( I ) ) 

YRC=PK*C0S I PHASERIT) ) 

YRS=PR*SIN(PHASER( n ) ^ 

Ync; = YMC + YRC 
YMS:iYMS + YRS 

' PHASER ( I ) “PHAS’rR'T f ) +OPHASERirj 

30 00 CO NTINUE 

3100 NC=:SIGMA*G.AUSS(0,0) 

^)S=SIG^;A*GAUSS(0«0 ) 

MC = rDC + Y^lC + MC 

ms=yds+yms+ns 

H( J)=SfiRf(MC**a+MS**2) 

DM< J)=SQRT{ YDC**2+YDS **g) 

*RH{0‘)=SQKT( YmC**2‘+Y>iS**2) 

CSO=SQRT( {Ydc+YHc) **2+(YdS+YHS)**2> 

IFtJ.EQ.JPU CSlsCSG 
IF< J.EO* JP2)CS2=CSG 
10000 COiJTIMUE 



c 

c 


* . * 

c ^ ^ 

C *_ calculate the response on the scan interval of the * 

C * candidate receiver ALGoRiTHM* * 

C * * 

c * * 


C *****************************************************^!i:4:^* 

C ■ ■ 

c 

C 

15000 THES1=THES 

CALL RCVRtMi JTtEST) 

C 

c ^ ™ 

c 

I eillGINAI/ PAGB JB 

THK(K)g TK OF PQQB QUALM 
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I 


i 


CS(K»=AMlintCSl,CS£) 

ERR-THETA-EST( JT> 

THER<K)=ERR 

ERMS=tRMS4-(THCR(Kn 

PRINT 1 b<Jo * EST I JT ) t THETA t THER ( K > 

15 pf, PORWAK ] Ho 1 / vl IXji * FJ.NAL ESTIMATES:*, El6 , S i sX » * TH£TAsit t £ 16 .fli SXi^ ERR 

10R-*iE16,ai’ 

S 

C 

■§ DECIDE IF S'h 6RT-TERM PLot“lS ‘OtSIKED 

C I F P LTSCAN IE TRU EiS .T. PL OT HaS ALRE ADY HEEN MaDE 

C 

_c 

IF (PLTSCAN) q.000 0 125000 

2SOOO C S K=CS( K ) _ __ 

CALL SHOUT (CSKiKi PLTSCAN) 

_C 

C IF SHORT RETURNS PLTSC AN=TRUE i A SHQRT-TERM PLOT IS TO BE MADE FOR THIS SCAN 
IF(PLTSCAl\n 26000, tfOOOO 

_c ; 

c 

_C ^ ORE ON FILE THE ARRAYS N EEDEO FO R THE SHORj-T E RM PLOT 

C ” 

26000 FLAG=,T. 

PRINT 26050 

26050 FORHAT<1HO,/,*A SHORT-TERM PERFORMANCE PLOT WILL SE MADE FOR THIS 
ISCAN*) 

E S T( JT+1)=THETA ^ 

□ o' 27000 1 = 1,5 

RA2( I) = RA(I ) 

RA(I)=RAi(l) 

27000 CO NTINUE 

TH£S2=THES 

thes=thesi 

CALL RCVR(DHtJTiESI) 

no 28000 1=1,5 

RA(I)=RA2(I) 

280 0 0 CONTINUE 

THES=THES2 

PRINT 261tl0,ESI ( JT) 

28100 FOKi'lAT(lHbtlOXi*ESTIMATE FOR DIRECT SIGNAL= *»E16.S) 

WRITE (20) FLAG 

WRITE(20) EST 

r.iRITE(20) M 

WR1TE(20) RM 

WRITE(20) ESI ^ 

URlTHaO) TA 

WRITE (20) DM 

60 TO 4'30 0 0 

40000 FLA6=iF, ^ 

WRITE(20 ) FLAG 

_C 

c 

JC 

C CHECK TO SEE IF THIS IS ONE OF THE 6 SCANS FqR WHICH A SHORT-TERM PLOT OF THE 

c received signal is desired 
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c 

c 

43000 IF(K.NE. 1NC8) GO TO 45000 

lFtK.GT,l+b*lNC) GO TO 45000 

C STORE ON FILE THE SAMPLES "OF THE 'RE'cEIVEU SIGNAL' CaNLESS THIS IS THf SAME 

C SCA M JiOR_l-JHlCH_THE_pn_A_lLEp SH0RT "TER H PLOT I S TO SE MAD E). 

iFl.fjOT. FLAG) WRITET20)' M 

IHCg=INC24-INC 

C 

C 


C UPDATE A/C and REFlECTOR COORDINATES FOR THE NEXT SCAN 
C 

c 

45000 KTK=KTK + \/*TUD ^ ^ " 

PRINT 1521iFLaG 

1521 fMTIAT (lH0'rroXt*FLAG=*iL2) 

theta=theta+dth*tud 

PHI=PHI+DPH*TUD 

IF< iHAX.Ea.O) 60 TO 50000 

do' 41000 I = lfIMAX 

RRTKU)=RRTKtI)+VR{I)*TUD 

THE TAR ( I ) sTHETAR ( I ) +DTH'Rri ) *TUD 

P4lR(I)=PHIRgHDPHK(l)*TUD 

4l000 CONTINUE 

sqooo continue 

ERNS=S0KI (EKMS/KMAX) 


STORTT^OfTFlLErn-HlTllKK^ NEEuTtrFOR THE lON<^-TERH PLOT 


C 

WRITE120 )ERMStTHK«CS.THER 
PRINT 1527 '■ 

1527 F0RMATciHi,9X , *SCAN K* i 12X t *THK ( K ) { MILLSEC M< t l3X * *CS f K ) * 1 17X » *THE 
1R(K) (UEG)*,/1 

DO GDOCO K-l.KMAX 

PRINT l323iK.TMK(K)iCS(K)f THERtKl 
l52S format ( llX.13i 3 (10X«E16.a) ) 

GOO'on CONTINUE ' ^ 

PRINT 1529, ERMS 

1529 FORMAT {1H6,//i 11X,*RMS ERROR- <^,E16.a) 

STOP 


END 


n n 


6 


i3!iLGAljS.S.i>LL 


: GAU SSIAM KANDOH wp . „GE11ER AJOK _ S .Hj., IHW _IN JRi 2/_25/ 75 


G "gauss utilizes the fortran RANOUH NO* function RANF{ ) TO generate RA'nUOn 

C NUMBERS UITH STANDARD GaUSSIaN UISTHISUT I ON ( ZERO MEaN. UNITY VARIANCE)* 


THE 1ST CALL TO GA USS IN T HE MAIN PR^^M SHOU LD B E WIT H X Ar AN IRRATIONAL 

C NUHBER(EG.tX=:sQRT(2.) ) TO INITIALIZE R'aNF. AlL OfHER CALLS TO GAUss 'aRE' MADE 

j:_wnii .x=aj 

c 

LOGICAL HAS 

iFtX.EQ.O,) GO To 1000 
A=RANF tX) 

HAS=. false. ■ ■ . . 

10 00 IF( HAS) 2000*3000 
2000 GAUSS=SAVEO 

HAS=. FALSE. 

RETURN 

3000 A=SQRT(-2.0*ALOG(RANF(0. n ) 

P=3. 141592705' *' ' ’ ' 

Tr2*P*^RANF(0, 1 

SAVEDsa'*SIN( i ( ^ ~ 

GAUSS=A*COS(T) 

HAS=',TRUE. ^ ^ " 

RETURN 

end 


I 







in n -n n 



S.H. .IRWjrjO, JH. . 7_/£5/75 


THft IS THE MLS SC AN AMGl £ C XM lTTLR ANTENNA 30RESI GHT ORlEM TATl OMt IN 

c degrees* where: 


C TS (SCAN TIME) =TImE IN SECONDS FROM BEGINNING OF tO SCAN UNTIL END OF 
C.^ SCAM 

C TR (KLFERENCE time") =TImE In S'ECONDS^ from WHEN THA =0 IN TO SCAn' UnTTL WHEN~ 

c^ THA:LD_DURING_FB0 scan. 

C TF rTlI^’E IN SECONDS FROM BEGINNING OF'fO UNTIL BEGINNING OF FWsCAN ~ 

C THo (THEtA-o) uVALUE of THA At beginning of tO scan 

c THs (THETA-s) =vaLue of THA at beginning of fro scan ' 

^ ^G sOMEG At t he scan RATE IN OEGREEs/SECOND. 

C ^ “ ' - - - -- 

c 


COMMON TR.TStTFiTHoiTHStOMGtTAU 
' T1=TS+TF 


IF(T.GE.O,0» go to 50 . 

THA=THo 



RETURN 



50 IF(T.GT.TS) GO TO lOO 



THA=THq+OMG+T 

RETURN 




100 IF(T.GE.TF) go to 200 


tha-ths 

RETURN 

200 IF(T.r7T.Tl) GO TO 250 
THA=THs-OHG*(T-TF) 




FlIMCTlnM PtTHETAl 


C 

C a nte nna selecti vity function 7/25/75 S.H* IRWlMtJR, _ 

C ’ 

C P(THETAi ISJTHE ANLENN^ SELECTIVITY AS A FUNgTIOM QF THETA i THE OIFFL'HEnCE 
C between SCAW angle and AIKCRaFT angular POSITION (A/C DEVIATION IN DEGREES 

C FROM SCANNING ANTENN A BQRE$IGHT), 

C 

C IF THE ~ABS'^Uix”VAXM~OF“(^^ ” 1 s“ L ES^ T HAN MI li. ' P'T A ) 'T^' S ^ 00“ P'l 1 *0 )sPI/4 

C (8Y L hosp itals RULE). 

C 

CONHQN/SLl/BEAHWTH 

MlN=l,OE~lO 

p I _1_A l5 g 2 6 5 3 6 

A=a,(|*fHETA'/DE;AMWTH 

B=1 . 0-A**2 

B1=A0S(B) 

IF(BI.GT.MIN) GO TO IDOO . 

P=PI/4.0 

R EXU_R N 

lOOn P'=CQS(PI*A/2*d )/B 

RE TURN 


END 


UBRQUTirJE yiHORTtgilUPS) 

LOGICAL BS 

dimension PiJ5) 

COMMON/Pl't/KMAX t ERR t PL 

— 

IF(K,£Q,17J PS=5,T, 

RETURN 

END 




1 


J 


I 


RUU KOUTI NE R CVHtMi J T^E STi 

50 UARE“G/\TE:. tracking receiver AUr70RlTHM S.H. IRWlNiJH. 9/10/7G 

REAL“n ' ' ' ™ ■ 

DINENSIQM W( JT ) «ES T1JT) 7RA{5) 

COMMON TRfTS.TF*THO»THStOHG»TAO - - 

COM MOM /TRACK/THE SiRA 

PI=3riAl5926536 ' ’ 

G=RA tl) 

del-0,0 

ALEST-0.0 



00 50Q0 Is:l»2 

JO-JP"10 

Jl=jp+10 

MN=:JO-H 

DO 2000 JsMNiJl 

DEL=0EL+0,5*TAU*G1*(H(J-1)*M( JJ ) 

EST( J)=THES+DEL 
IF(J.NE.JP> GO TO 2000 
61=-G1 

ALE ST =ALEST4-M( JPl 

2000 CONTINUE 

PRINT 105, JOi JPi Jl,DEL - 

105 F0RMAT(l5x»*J0 = *t I5,5X»*JP=*»I5t5Xt*Jl=*»l5.5X»*DEL==t‘*E16.6) 
IFtI.E0,2> GO TO 3000 
DO 2100 OsltJO 

EST(U)=THES1 

2100 CONTINUE ' ' 

MINsJl+1 

JK-’Jl 

5000 JP-IFIXC tTF+(THS-THES>/0MG)/TAU+l,5) 

Gl=-G 

5000 CONTINUE 

THES=EST(dl) 

D EN=COStl,2*PI*OEL) 

lF(ABS(DEN),LT.1.0E-e) GO TO 5100 
ALEST= ALEST* « 1 . 0- ( 2.4*DEL ) / ( 2, Q*DEN ) 

nACns-ONG/ {4.0*ALESTT“ " 

5100 DO 5500 JsMiNfJO 

ESt ■( JT=EST ( JK ) “■ 

5500 CONTINUE 

MN=Jl+l 

DO feOOO UsWNfJT 



ESTlJ)=EST(Jai 

6000 corm.JuE 


RETURN 

END 



n ooln nlono I o n o n r».n ri p o n n 


_^R_0 G P /\ M L,S P^L pT { I N P U T j.Q UJ PU TtlAP ■ TAPE13 . TAPE151 


T/VPE12 = PL0T PlL't TAPEl3=PERMAf4ENT PIlE 

JHIG PKOC’I'Af-i PL-PTS THE RESULTS OF THE pROGRAn nlr.SRpXR i >'HTC}i ,ARE 

stored FOK a short time on PEnMAPiENT FILE. ALL INFORMATION NEEDLO BV 
_JH 1 5 _ PJ^ 0Ct« AM T S_ TQ. BE „REAp._FjRO,M THE FI LE I N UN FOR MATTED FORM. THE 
ORDER Is" AH FOLLOWS. A SLAsH DENOTING ENO-OF-RECORD : 

JOPTjiNAME,NRCVRjCHmN,ERsPECt JB.DAY/ 

kma X 1 1 MAX , Alpha , alpham . snr i sigma . ti » jt/ 

RTK, VjJHETA.DTH.PHI iDPH/_ 

RRTK(l) . VR (1 ) .THETARO. > .DTHR (1) . PHIR ( 1 ) . DPhR ( 1 ) . ALPHAR(I)/ 

RPjKtS) , VI; ( 2 )',iHEj[ ^RL2J .0 tHR( 2) . P HIR ( g ) . DPHR t S > . A LPHAR (2)/ 

“rRTkT3') , VP ( 5 ) 1 tPE'. AR < 3 ) , DtHR ( 3 > , PHIR ( 3 J , O'PHR ( 3 ) * ALPHAR ( 3 J / 

TM/ jscAw Time array) _ 

■■'o AT a" for' detailed 'and "sFLECTED short-term plots/ (EACH SCAN IS 

REPRESENT ') BY A FLAG. iF FLA^= .TRUE._t_A_ DETAILEJ) SHORT TER M PLO T 

*’lS TO BE MADE. PLOTS OF THE COMPOSITE SIGNAL WILL BE M’AdE FOR SIX 

uniformly spaced SCANS. BEGI NNING WITH T-HE FtRST.) 

ESTERR.TMK.CS.THER/ 


_ logical FlAG^ 

REAL NOISE .Ml. M2 

_UI MENS ] orj ,cs (^0 )j THER ( 30) . T MK(30) .NQISE(3o) . TM( gOOO ) . Y ( 20 00 ) 

dimension ERRPCA<4) .A> ^0S(5) .SR1P0S(5) .SR2P0S(5) tSR3P0S(5) 
dimension ACRATE(5) .SIURAT(‘f) .SRgRAT(4) » SK3RAT ( 4 ) » DP AMPA ( 2 ) 
"dimension RlAMPA(2j .R2AMPA ( 2i , R3AMP'A( 2') , SmA (> ) , sNRDBA (X)", KTRA ( 2 ) 
OIJIENSlOlj ESTERA( 4 ) .NRCVR (3 ) . KMAX AJ 5J i„STLbL ( 3J . TRANGL ( 3 ) . NAME ( 3 ) 
"OIMENSICN RRTK(3*) . VR ( 3 ) VTHETAR ( 3 ) . DTHR ( 3 ) .PHIRf 3 ) »DPHR (3) iAUPH'AR(3 

1 ; .NRCVR2( 3) 

CALL CALC0MP(i2) 


READ INITIALI ZATIO N WTA .fLR0H._?lLE 

REA0( 13) I OPT. NAME. NRCVR. NRC VR2.C HAN . ERsPEC.JBv PAY 

READ ( 13) ‘kMAX, iMA)f.' ALPHA.'lALPHAH.SNRisiGMA."Tl. JT 

READ (13) RTK.VjJHETA.DTH.PHItDPH 

■ ■■ DO in 1 ^ 1 , 3 ' 

RE AD { 1 3 ) RR TK { I UVR < I ) , THETAR ( 1 ) . DT HR ( I ) «P H IR ( I ) .OPHR( I ) . ALPHAR ( I ) 
ih CONTINUE 



i 


i 

- 1 


J 

n 

■ i 

J 

■ I 


ENCpDE_P_LriT labels ; j 


ENCqOE(27.5.KMAXA)KMAX 

ENCOOF ( 1 R , 1 Q . DPAMPA ) ALPHA 
ENCODE lIO . RIAMPA ) ALPHAR ( 11 
encode ( in , 10 . P2AMPA ) ALPHAR I 2 ) 
encode n (■ , 10 . K3AMPA ) ALPHAR 3_1 
ENCdbEl lAtl O.SHA) alpham" 

ENCODE ( , 15, ACPOS ) RTK .THETA .R4J 

ENCOOU4 2,2 0,'aCRAT£)V'.DTH.DPH 

ENCODE ( ‘t 1 , 25 . SRK’OS ) RRTK.(_1 ). 1 1.HLTAJ (1) .PHlR(l ) 

encode (37,30. S)UHAT ) VR ( 1 ) .~DTHR ( 1 ) . DPHR fl ) 




E:}\JC 0 DE: • 25 _t_SK 2 POS)RRTfv( 2 ) t THF_TAR( 2 ) iPHlRig) 


] 


o n o n n n o nio n n ^ ! In n'n r> ' ^ on n 


I 


encode: < 37»30, SR2RAT ) VRt?) iDTHH (2) t OPHR t 2 J 

FNCODE:(41_*2S»SR3POS)RRTK(3) .THtTA^( 5) *PHl R(3) 

ENC6DE(37,30tSR3RAT) VR{3) t OTHR { 3 ) , DPHR { 3 ) 

SN2=SIGMA*SQRT(2.0) 

encode (6f35iSNA)vSN2 

__ ENCOuE(16.40*SWRDBA)SNR 

encode ( 40 1 50 I ERSPCA) ERSPEC 

5 FORMAT (*U0NG-TERM .PLOT FOR SC ANSt) 

10 FORMAT r*AHPLITUDt: = *,'F6.2) 

_15_F0RMAT ( tl-VITlAL- A/C POS_ITION = *tF6',2t*« *_tF6.2i** 

20 FORMaT(*INITIAL A/C RATE = *tF7»2i*» *tF6.2»*t ■mF6,2) 

25 FQRMAT(*INITIAL POSITION = +*F6.2t*t * * F6.2j * » _ 1 1 Fg, 2) 

30 FORMAT (^INITIAL rate F6. 2, * t *i F6 .2 , + , . F6. 2 ) 

35_F_0KMAT (F6.2) „ 

40 FOHMA'TUONR =5 *iF5.2i* DB)*) 

50 FORMaTURMS error SPECIFICATION ^ *> F6.3<» pEGREES^t^) 

“read “scan flME^AXlS "and" scale FOR* PLOTTING *“" 

^dTX3) TM 

CALL SCALE(TK,25.0. JT,1) 

VMIN = TM((JT + l)*100b’ ' 

J1£L T A = Tj;i{ J T + 2 ).* 1 0 00 

I N I T i*AriTr“^FsELECTE^ sHORT-TERM PLOTS 
TncskmAxTS — 

__ KPLT = 1 

“*D0 GO 0 0 KsltKHAX 

R£AP< 1 3) FL AG 

1F( .NOT. FLAG) GO" to 2000 

KTRUe=K __ 

IF(K,Ne.KPLT) go to" 4050 

kj>lt=kplt+inc 

4050 IF (10Pt J'J*E^2') *G*d TO -5000 

^ 4500 KS = 1.6 

REAU(13) Y 

4_50_0 cONTlj^UE 

GO TO 60 bo 

2000 IFIK.EO.KPLT) GO TO 3DD0 

' go' to'gooo' 

3000 IF( K.GT.5<=INC4-li G O TO 60 00 

RPLT=KPLf+INC 

read (13) Y 

GO “to 600 0 


_♦ ♦ 

* * 

* GENERATE DETAILED sHORt-TEKM PLOT * 


****** 4^* + **+**Hl * + ♦;*;** lt*#**5t!* + *****.*!**^;*^!* 




c 

c PPj9JL_Ea££5.J^IlQ^AI^piNG . harks 

C 

^Q.Q.Q-..L5k5l 

GO TO G2000 

6£10.0 -CJQmWL^ 

REflDdS) Y 

mLlAX=Jf-< JT+11 

CALL SCALclYia.O* JTtl) 

ELSXJ.:=y_( JT, + 1 ) ^ : 

EST2=Y( JT+2) 

aEAIl{13) Y : 

call SCALEI(Y»2.0t JT»1 ) 

S£1£Y1JT -ai 

READ(13) i 

^CALL_SCALE( Y*2.^tJj,l >■ 

S22=Y(JT+2) 

READ<13) Y _ _ 

CALL SCALE(Y.2.6t JTii ) 

e:s ii =yijtuJ 

E;SI2=Y{J7+2) 

RMOU3) Y 

READ<13)Y 

q A LL SC A L E i Yjl2 .OiJT . 1_L. 

S23=Y( JT+2) 

Rl = 0 . 0 

M2=AMAX1(S21»S£2»S23)*2,0 

BACKSPACE 1 3 

'"‘backspace 13 
BACKSPACE 13 
“CACKSPACE 13 

BACKSPAC E_J. 3 

BACKSPAC'E 13 
KEAD(13) Y 

IF<ESI2,GT.eST2) GO TO 5110 

ESIl^ESTl 

ESI2=EST2 

GO TO 5 1 30 

5110 EST1=ESU 

E ST2= ESI2 ; 

'EST2=ESI2 

_5130 Y_f JT+1J=EST1 

Y(jt+2)=EST2 

- C - 

C 

C D R AW THE R ECEI VER OUTPUT FROM THE C OWPOSTTE S IGNAL _ 
C ' 

C A L L_ R, 0 ? { 0 . 0 t0t2»“3) 

'S = -ESTl/EST2' 

I£i ESil.. GJ , 0.0 ) $=Q . O 

XMAX=Tl/Tin( jf+2» 

XMl-XMAX+0,1 
XK2=iXMAX + 0.'2 ■ 

C Aj-^LOT ( 0 , 0 t ^3J 

CALL P'LOT{XMAX*»St2)' 

CALL AXIS(O.OvO,OilH 1 1 » 2 « 0 i SO . 0 ♦ ESTl t ESTP ) 

SX=(Y<1)-EST1)/EST2+0,1 

C ALL S YMOOL{Q.ltSlt D. lAt 1A H( DEG, ) i 0.QU4 ) 


IF(Y(1),GE,1,0) GO TO 5134 
CALL SYMDoL(0.5iSlt0,14.i lHO»O.Otl) 

5134 C AUX NUMBER ( 0 . SlVo , li^Vv ( 1 M 0 .'o t 4HF7 . 2 » 

JL_ESI,. 

CALL Llh}E(TM»Y»JT*i»0»0) 

S15{ Y (JT)-ESTl I/EST2 _ _ 

CALL SYMBOL<Xrn*Sli D.14a4H( "‘DEG", ) , q", qV iX) 

IF t Y ( J_TJ, .J5 E . 1 . 0 ) GO^J)„ 513 5 

XM3=XM2+0*5 

call SYMBOL(XM3tSlt0.14flHOtO,0»l) 

5135 CALL NUMBER (XfTgt Sit oVl't-* Y (JT ) .0,0 t 4HF'7. 2) 

CALL SYMBOL(-2.5,l,56t0.14t9HCANDJOATEto.0.t?)_ _ 

CALL SYMBoL(-2*5,i.26t0.14tllHRCVR OUTP JT , o'. 6 1 11 j 

call SYMBoL<-2«5.1.00t0.14t9Ht0EGREES) to.otgi 

CALL SXMadL<-2.5tO,5&i 0.14il6HCbMPOSlTE ‘SlGNALtO.Oi 16) 

CALL SYMBOH-2.51 0.30' 0»l4tl0HPLUS NOISE* 0 . 0 * 10 ) 

c . - 

c 

C draw the COMPbSiTE SIGNAL AMPuitUOE VS TIME PLOT 

. 

rEAD(13) Y 
Y( JT+1)=M1 

Y(JT + 2)'=H2‘ ■ * 

CALL PLOT ( 0. 0* 2* 5*-3) 

CALL PLOT UMAX '0.0 i2) 

CAL L AXI5(0«0 , 0. 0, I H tlt l.O , 90 , Q , Y ( JT+I ) , Y { JT-fg ) ) 

C H 

CALL LiNE(TM.Y,JT,1»OiO) 

CALL SYNBoL(- 2.5*0.7* 0.i4ii2HrA'HPLlTliOE OF,b*b,lg) ^ 

CALL SyMB0L<-2.5,0,42, 0.14,16HCOMPoS1TE SIGNAL ,0 , 0 , 16 ) 
call SYmBOL(-2.5,0.14,O.14,10HPLUS n 6IS£ ,b; 0 , 10 ) 

C 

C 

C DRAW THE MUlTIPATH AMPLITUDE VS TIME PLOT 

'C ■ _ : 

READ{13) Y 

Y'(JT + l)=Ml ■ ' _ 

Y(JT+2)=M2 

CALL PLoT ( 0.0,1 »b«-3) 

CALL PLOT(XMAX,Q.O ,2) 

CALL AXIS(0.0 .O.OVlH , 1 , l". O', 90 , 0 , Y ( jT+i ) , Y t JT + 2) ] ■ 

__j: RM 

CALL tliNE{TM,Y, Jt*l,0,O) " . ' 

CALL SYMB0L(-2 ♦ 5 , 0„.,lLgj, l4 , 12HA MPLIT UDE OF, 0,0,12) 

CALL SyMB 0L(-2. s', 0 . 42 X 0 ^ 14 , 9HMULTIPATH, 0,0, 9) 

CALL SYMB0L(-2.5, 0.14,0.14,12HXNTERFERLNCE,D.O, 12) 

“c " " ■ 

c 

C DRAW THE RECEIVER OUTPUT VS tlME PLOT (ONLY DIRECT SIGNAL PRESENT) 
C 

— CA UL PLOT( 0» 0 , 1 ^ 5 , “3) - - 

READ(13) Y 
y(JT + l")=ESir 
Y(JT+2)=ESI2 
S=-ESI1/ESI2 

1F(ESI1. GT, 0 , p )_J=0 , 0 

CALL PL0t(d.b,S,3) 

CALL PL0T(XMAX»S,2) 


1 


CALL AXIS{D.0*0,0tlH il * 2 . 0 1 90 . 0 , ESIl *esi2 ) 

SJ -J IUJ_C.E^I 1 1 /E.S 1 2+il , 1 

CALL SYMB0Lt0,ltSliD,14,14H( DEG * ) » 0 » 0 1 14 ) 

JFiY(l).GE,l.n> GO^TO.5138 

CALL SYMDOL(0.5tSlfO.i4ilH0id,0il) 

_&13a_CALL WU(^B.ERjO.£vS 1»E...14*Y(1) »0^0.l.4}^F7.2^ 

C ESI 

C ALL_LJNf;iTM » Y» JTtltmn ) 

S1s(Y(JT)-ESI1}/ESI2 

CALL SYMBOL (XHl *_S1» 0 , 14i 14Ht D EG , ) , 0 , Di 1 4 ) 

IF(Y( JT) .GE.l.OJGO TO 5136 

X M 3s X M 2 + 0 , 3_ . 

CALL SYnBOL{XW3«Slf0.l4,iH0iO,0tl) 

5136 C ALL NUFIBER( XM2 , ,S_^i 0_.Jj+ , Y ( JTJ. , 0 ,0,4H F7 . 2J 

CALL S y'mBOL (-2,5,1.56,0.14 ,'9HCaND IDATE'i 0.0,9) 

CALL SYMBOL(-2,5,1.2R,0.14,iaHRCVR OUTPUT , 0 . 0 , 11 ) 

CALL SYMD6L(-2»5,1.00, 6.14»9H(DEGREES) ,'o.Ot9Y ' 

CALL SYMB0L(~2.5, 0.58, 0.14, IIHDIRECT PATH, D. 0,11) 

CALL SYMB0L(-2. 5, 0.30, 0,14, 14HCOmP6n£WT Only, 0,o»l'4) 

X ^ 

c 

C DRAW THE SCAN ANGLE AND DIRECT SIGNAL AMPLITUDE VS TIME PLOTS 

C ' ' 

CALL PLOT (0.0, 2. 5, -3) 

READ(13) Y ’ : 

CALL SCALE(Y,2.0, JT,1) 

S=-Y('JT + 1)/Y( Jt+2) _____ 

CALL PLOT(0.0,S, 3) 

CALL PLOTY XMAX » S , 2 ) 

^C^L_AxIS(0,0,0,0,1H ,1,2.0,90.0, Y{UT + 1) , Y(JT+2) ) 

C TA ~ ' "■ ■ “■■■ 

CALL DASHLNCT M,Y,JT, 1,0,0) 

CALL SYMBoL{-2'.”5,i.5,0.14,10HS.'.'AN ANGLE , Q”76", 10 ) “ 

call SYMB0L(-2. 5, 1.02,0. 14, 9H(DEGREES) ,0.0*9) 

CALL symbol i -2 , 5 , 0“. 6 ,0 .14, SH ( DASHED ) , 0,0 ,T) 

CALL PLOT(0.0,S*-3) 

READ(13) Y 

Y(JT- )- 1)=M1 

Y(UT+2)=M2 

C DM 

CALL LINE(TM,Y,JT, 1,0,0) ^ “ 

CALL AXIS UMAX, 0.0, IH , -1 * 1 . 0, 90 . 0 , Y ( JT + 1 ) ,Y( JT+2) ) 

CALL SYMB0L(25.25, .65, 0. 14,12HAMPLffUDE OF, 0.0, 12) 

CA LL SYMBOL ( 25 . 2 5 3 7 j,0 . 14 , llHD J RfCT PA TH, 0.0, 11) 

CALL SYMBoL( 2 5.2 5,0, 0"6,0.l4,9HCOMFQNENf, 0.0,9) 

CALL SYMB0L(25.25,-.25, .14,12H(S0UID LINE ) , 0 . 0 , 12 ) 

c ^ 

C 

C SET ORIGIN an INCH PAST LAST FOLD ANO AT Sotf d«" 0F“PAGE 
C 

ADJUST=8.9+S 

CaL L PL OTT 26 . 5 , - a D JU S T , 3 > 

C ■ '™~ 

c LABEL Tiftir^c'^K ^ ORIGINAIi 

E N c 0 dE(27,510oTstlbl)k OF rOOK 

5100 fqrmat ( *short -ter m Plot f or sc an *,i2) 


I 


1 


1 


KTK=(K-l)+75 

ENCODE ( 15 ,515 0 t STLBL2 J KTK 

5150 format <iTK :: *,m,*MSEC. + j 

ENCODE (29, 520 0iTRANGL)TMETA 1 
5200 FORMAT (*A/C AT **F6,2«* DEGREES *A2 iMUTMi'r’ 

_ CALL SyMOOLd, 0,10. 22,0, 14, STUDL,0,0,27}__ 

GO TO 64000 

6 410 0 CONTINUE ^ . 

CALL symbol ( ,75, .665, 0.14 iTKA'nCL, 0,0,29) 
CALL SYM0OL(O.75,O,586, 0.14,ST_LJL2j0,0ii§i 

6000 continue 

IF(IOPt.EQ,3) go to 17000 

READU3) ESTErK*TMK,CS,7HER 

I FJ 1 01^ EQ , 2 ) _G0_T'^IM.(IP 

CALr'MEWBLOK 

BACKSPACE 13 

" KPLT=5»INC+1 


**** :,:**** *******:iK ************* t*******;****^: ** + *:!« *!*#***!!:***:***** 

* * 

* * 

* G E N^_A_T E _S E U EC TED SH OR T‘»TE R|^ PLOT S pF THE RECEIVED SiG N A L * 

* * 

_* * 

***************************************************************** 


LBL=2 : 

Nl-1 

_G 0 T o 6 2000 

continue 

call PLoXfO.O, D,2t-31 

"bo' 15 00 0 KC 0M= 1 , KH AX 

_K5KMAX+1-KC0M 

FLA'Gs, false. 

IF(K. E S.K TRUE) GO TO 6000 ; : 

IFCK.NE.KPLT) g"0 TO 14000 

KP L T = KP L Tji^I N C 

BACKSPACE 13 

_P_EAP(13J_X„ 

IFiK.EQ,5*INC+;) GO TO 12500 
CALL PLOT{Q,0, 1.75,-3) 


CALL PLOT(25.0,0,0,2) 

IF(NJ^.NE.I) GO_TO_12510 

"call SCALE( Y,i.o,'jt 7 u 

AM 1 - 0.0 

AM2=Y( JT+2J 
Nl=2 


y (JT+l)=AMl 

jrj JX+ 2 ) = A M 2 

CALL AXIS(0.0*,0.b,fH , 1 , 1 . 0,90 , 0 » Y ( jt+1 ) , y ( JT+2 )?' 

CALL LINE(T«,Y, JT,1,0_,U 

"backspace "13 

IF{ .NOT . FLAG) G O TO 14OC 0 

BACKSPACE 15 

backsp ace 13 


15000 

c 

CONTINUE 

c 

..C .LABEL AXE.S 

C 

DO. 65&50 NK = 1,6. 

NSCAN=1+(NK-1 )*INC 
.ENC_OD.E(7vS5010»NSCANA)NSCAN 

65010 

FORMaT{»SCAN ♦iIB) 
KTK=75*(NSCAN-1) 


ENCODE ( 12 1 6502 0»k‘TKA)KTK 
. format {*Tk = *,I4,* MS*) 
call SYMBoU-2.5.0*57»0,i4iNSCANAibt0t7) 
CALL SYMOQL(-2.5iO,29tO,l4*KTKAi 0 ,0tl2) 

IF(NK.E0.6> GO TO 65050 
CALL PLOT(0.0,-1.75,-3) 

65050 

CONTINUE 

CALL PLOT (2G.5t-0.9,-3) 

CALL SYMBOL(fl,0,10,22,0.14,40HAMPLITUDE Op COMPOSITE SIGNAL PLUS N 
lOISE.t.O.OmO) ... 

64200 

GO TO 64000 

continue 

CALL NEWBLOK 
GO TO 16000 

6000 

IF(K.NE.KPLT) go TO 9000 
flags* TRUE-, 

KPLT=KPLT"INC 
DO llOOO kB=1,5 

liOOO 

BACKSPACE 13 

continue 

9000 

go TO 12000 
DO 10000 KB=1,6 

10000 

backspace 13 
CONTINUE 

c 

GO TO 14000 

c 

c 


c 

c 

* * 

C 

c 

* * 
* generate LONG-TERM PLOT * 

c 

c 

* * 

* * 

c 

c 

« « * :|^ « * * * * If: * * « « * « # 4: * * « * * * 4: 

C 

16000 

LBL = 3 

c 

C FILL NOISE ARRAY FOR PLOT 

C 

DO 66000 KL=1^KMAX 

66000 

N0ISE(KL)sSrGMA*SQRT(2,D) 

CONTINUE 

C 

c 



C DRAW THE EDGES AND FOLDING MARKS 
C 


A-22 


CALL PLOT(0,0,-0.5»3) 

CALL PLOT(O.OilO,bi2)^_ .. . 

CALL SYKnoLfa.bi-.HtOVait 13i D,0i-1> 

CALL SYhtJOLCft, 5*10. 4,0,21,13.0.01-1) 

call SYHB0L(14. 5*-. 4, 0.21,13,0, 0,-1) 

CALL SYMBOLll4.5,l0,4,0,2l,l3,0.t0t-IJ 

CALL PLOT (20. 5,-0, 5, 3) 

CAL L P LOT (2 0 . 5 , 1 0.. .5., 

C ' ■ • 

C _ 

C move origin to 3/4 INCH FROM CUTTOM OF PA>ER 

C „ 

c 

CALL PLOT ( 0, 0.25, "3) 

C 

^labe:l_top axis 

c " 

call SYMBOL(1.5»6.35,0.14,16HAV. AMPUITJUe; OF_,_0 , 0 , 16 )_ 

'call SYr'iBOLd, 5, 0,07,0.14, 16HC0HpbsiTE SlgNaL , 0 , 0 » 16) 

call . S Y Fi H 0 L ( 1 , 5 , 7 , 7 9_,.P. tl 4 , .1 7 H A,T_ T. I E Q F_ DiB L C T ,i_0, , 0 , 1 71 

call SYNBdLll,5,7.51,0.1‘),15HSIGNAL CEWtROID, 0, 0 , IS ) 

C;LL SyMBOL(1,5,7.£3j0.14,1£H(S0LID LINE) ) 

'■-.’■^L SYMB0Lll.5,6V67,0‘.i4,15HRr>;s No'iSE uEvEL, 0 . 0 , 15 ) 

C ALL _.S YMBOL( 1 , 5 ,6. 3_^1, 14, 1 3H(DAS HED LIN E ) , 0 . 0 , 13) 

C 

c label b o ttow axis 

c 

_CALL SYMB0L{1,_5,2,1, 0j^4, IGHESTIN^TION ERrOR , Oj^O » l6 > 

'CALL SyMS61V1.5V1.‘62,0,14,16H(A/C' ANGLE-Est. )‘‘,0.b,l6) 

CALL SYMBOLll . a,l • 54, 0.14 , 9 H_<^e:6REES_)j 0 ,0,9) 

c 

c nQVE origin an Inch to right o f lettering and up 

c 

CALL^PLOmj 5,5. 3J 

C ' " 

c scale t he data for the error plot 
c 

CAL L SCALE( TNK,7.0, KPIAX , 1 ) 

CStKMAX+Dsd.d 
CALL SCALE{CS,4,0,KMAX4-1,1> 

CSCKNAx+'1)sCS(KMAX+2)' 

CAtKf'lAy + 2)=CS(KMAX+3) 

NOISE(KHAX + l )=CS(Kha'xVi) 

NOIS E(KMAX -t-2)=C S (KHAX4-2) 

C 

C DRA N AXES ^ 

C 

CALL AX1S(0.0,0,0,35HTIP'1E(SECL SA^CE S TART OF FI RS T SCA N, -3&, 7 . 0 , 

lb.Q,t«K(KMrX+l) ,TMK(K«AX+2) » 

CALL AXIS(0. Q, 0.0 , IH ,1, 4 . 0 , 90 , Q , CS ( KMA X-1-1 ) i CS ( KN A X-l-2) ) 

_c 

C DRAW PLOT WITH AM ASTERISK AT SELECTED POINTS 
C 

CALL LINE{T^!K,CS,KNAX, 1,0,0) 

If ( IOPT.EQ.2) goto 3005 

IMX"5*INC + 1 

DO 3001 "isf, IMXTiNc 

SX=(TMK{I)-THK(KNAX4-1) 1/T»K(K HAX + 2) 




u-u 0| ■ ookj k>u o I «_* 


S1 = (CB(I l-CS(KMAX + l) )/CS(KI^AX + 2) 

__C.ALL„SYMBgL {SX i Si t 0,, U 1*0. 

3001 co.niNuE 

SX=(Tl^K(KTRUr )-TMK(KMAX+1 > ) /TMK ( kMAX+2 ) 

S1=(CS(KTRUE) -CS(KMAX+1) )/CS(KHAX+2) 

CALL ..SYMBOL tSX . Si . 0 ,.ltf » 1 1 0 , 0 1 -J ). 

3005 CALL DASHLN(TMK«NQiSE:iKMAX» ni)*0) 

c 

C MOVE ORIGIN BACK DOWN 

CALL PL6f(0.0*-5t5*-3) 

c” Y calT Error” bYta "" 

CALL scale rfHER riT^o"! KMAXTf) 
cYYf^TC”TR'uFYRISlN " 


TOR=0.0 

LFeTHERLKl^X+l L.GEj_QJ GO TO 5111 

T0R=-THER (KMAX+l ) /THER ( KMAX+2 » 

3111 C O NTINUE 

DRAW A XES 

CA LL AXIStO.f i » TOR t SSHT iME ( SEC ) SINCE START OF FIRST SCAN»" 55t 

17.0i 0,0»TMK (KMAX + l) , TMk‘(KMaX + 2) ) 

CALL AXIS ( 0 .OtO .Oi iH .D»9 0»0.T HER ( K*1AX4-1 ) * THER ( KMAX-i-2 > > 


DRAlj PLOT W IT H AN A STERI sK AT SE LECTED POINTS 


CALL LIN£(TMK»THER,KMAXtli0t0) 

IFt lOPT.EQ.2)' Gd TO 3006’ 

DO 3000 I=ltlMX»lNC 

S’X=(TMK ( I )-TMK IKMaX + 1 ) )/TMK (KMfiX + 2) ™ " 

S1=<THER(I)-THER(KMAX+1) )/THER(KMAX+£) 

CAXU SYMBOLrSX»SliQ.lti»ll’*’d.O*-l) 

3009 CONTINUE 

SX=tTMK(KTRUE)-TMKtKMAX+l) ) /TMK ( KMAX+2 ) 

S1=(THERIKTRUE)"THER (KMAX+l) ) /THER < KMAXf 2 ) 

CALL SYMBOL (SXrs’ltO.mt It 0.0'»-1) ™ “ 

move origin and T r'e ^e “t iTl’e ' b l oc k ' 

3006 CALL PLOT{9.75t-0,75t-5) ~ ~ — — 

C A.L.L_1YMB„0 L(1 . 0 1 10 . 2 2 t 0 . 1=^ . KMAXA , 0 . 0 1 27 ) 

GO TO 64000 'r 

64300 CONTINUE 

encode ("36 1 45 ^ e'sYe’KA ) ESTEITr 
45 FORMAT ( *RMS ESTIMATION ERROR = »,F5,3i*: DEGREES*) 

CALL SYMBOU ,75t .665', lAi’ie'HESTlMATIO'N ERROR t o”, 0 * 16 ) 

CiU.L_ SYMBOL t .gSt .655. O.lOStESTEHAtO »0t 36) 

CALL ’s’YMBOL r. 95 1'. 4 45t 0 . 1 0 5 , ERSPC A i 0 3 , 4 6 ) 

IF( I0PT.EQ.2) GO TO i7000 

CALL SYMBOLt ,75. .i65t 0.07Vxit ,0i-l) 

CLA.-L_SXM80L ( , 85 t_._lG5ji 0 . 0 7 1 4 SHSH ORT-TERM Pi QT OF COMPOSITE SIGNAL A 

lVAiLABLE»0.0»45> 

CALL SYMBOL (0.75* ^,_0 2 5^ 0 . 0 7 1 O.Ot-1) 


I 


CALL S-"!iOHP.05*0.a?5*0.07*4..mjt:TAILE:D SHORT-TLRM PERFORf^ANCE. PLO 

IT AVAUAUlE < 0.0 t'K, > 

GO TO 17000 ' ' " 

C .. . . 

c 

c . , 


j; * *. 

C ^ ' 4t 

_C _ * DRAW THE_ EOGEsS AND FOLDING MARKS * 

C * * 

c * _ * 


c * + ***>♦ i»***+*>{=*t**+*j|>****<t*#**!»:**:tl*lM:*** 

_C 

c 

c ^ _ 

6?Oo6 call PLOT(0.Gt-0,5t 3) 

CAUL PLOT ( 0,0 . 10.5» 2 ) _ 

CALL .SYMBOL(C,5t“O.Ai 0.2itl3t0.0.-l) 

CALL SY''‘t^OL{fl.5.10.4,0.21,_13.0.0i..ilJ 

CALL SY-:bOLtl^.5,10.K,0,21,12t0.o'*-l) 
call 5Yf'HiOL(15.5.-0.‘+t0.21tl3t0.Qt-l) 

CALL SYMPOL{22,5'~0.l. 1 ,21.l3»0.Ct-l) 

CALL SYHiOL{?2.5'L0.K*0*2l,l3i0.0t-l) 

"' call SYi':P.f)U2'?.b,10,4,t).21,13,0,oV-lj 

__ _ CALL„SYF(.‘3Lt29,5i-0,4.0j,Jl_il3»O.G*-_lJ 

“call plot (36. 5. 10. 5^3 > 

call Plot ( 56.bt-0.5<2) 

"c 

c 

c tmAl,) the tike axis 

c _ 

CA'lL plot (i+.b.0.2*-3) 

call axis(o.o,o,o,3?hreceiver Scan timE(mtuliseconps) ,^ 32 , 25 . 0 , 

'''lO.'OtV';iNt;)ELTA) 

G(^T0( 62100 .G?200)LliL __ _ 

c'“" 



C * + ***4:«.**+;*i;4 + 3ti*:^ + **4H;**^*****:!(:** 


c __ *__ __ ^ 

c * ^ 

C *_ jlTLL BLOCK BODY FOR PLOTS 

c ' * ' ' *’ 

C t » 


c 

C ■■■ ' ■ ■ "■ 

64000 CALL SYHBOUl.Liy.OL^ 0 , 14 i 5 HDATE: tO.OtS) 

' CALL SY:-POLt2. 11, y. 94,0.14 , DAYtO.O.'IO) 

^ALL SY iLOL( 1 .4 , 9.h6, 0. 14 , IbHPKOGRAM! HuSRCVR , 0 . OU6 ) 

CALL SYi'-|‘OL(1.4,9.3a,0.14,9HJOHNAME; iO,o',g)' 

CALL SY'-BOL(P,43,9,38.0.14,JB«O.Oi7) 

'CALL" FYKr, 0 L( 2 . 10 , 9 . 10 , 0 . l 4 *NAMEt 0 ^bi 3 O'r'' 

CALL FYTPOL(2.10,9.1fj,0.14,NAMt ,0.0*10) 

CALL SYYBOU .75,8.541 0. 14, QHCHANNELI* 0.0,6) 

_ CAL L SYKMOL {1.93,6.54,0.14, CHAN ,0.0,10) 

CALL“SYHIiOLV,75,6;12,0‘.14',I4HS1cNAL MAKEUP ,0.0’,14) 

____ _ _C^LL SYMliOL( . 75 , 7 . 96 , 0 . 07 t 47 HPOSITION SPECIFIED BY RANGE I N NAUt IC 


A-25 


lAL MlLESt t0.ntA7) 

CALL SYWB0 L( . 7bi 7,84 » 0 . 07i 41HpLGRECS,.pF_A£IHUTHj ^DEGRrCS. pF^ELEVAT 

1 1 Oh\l , 1 0 . 0 t 4 1 ) 

CALL SYMHOL( .7bi7.7t0.07.46HRATES SPECIFIED BY V/ELOCITY 1 KNOTS* 

1 OEGHEES ' *0.0*46) 

CALL._StnUpLt ,7b *7,56* 0.07* 35 HAZIMUTH/SECi„DEGRELS ELEVAT J ON/StC , , 

10,0*35) 

C ADJ UST ORIGIN . 

CALL PLOT (0,0*0,5*-3) 

CALL_SYMBOL(. 95*6, 78,0, 14*11H0IRLCT PATH , Q., 0 , 1 j, ) 

CALL SYMBOLa, 05*6, 57*0. 105*DPAMPA, 0.0*18) 

CALL SyWIBOL<1,05,6,36,0.105*ACPOS*0.0*45) 

” ’ CALL'SYMGOLn,05,6,15.0.105,ACRATE*0,0*4l) ’ 

ilALL SYWBQL( .95*5.765* 0,14* 16HSPECULAR 'I UlT I F ATh * 0 , 0 * Ifl ) 

lV(iMAX.*GE.l)GO TO''64010 

C ALL _SYMB0L (3.2 5 , 765 * 0^. 14 * 5H-N0NE , o 5JL 

GO TO 64020 

64 010 call SYFIBqU .95*5. 465.0. 1225.11HREFLECT0R 1*0.0*11) 

CALL SYMBQL(l,05*5.275,0.10S*RlAMPAtO.0,l6) 

C ALt^SYMBOLJl. 05 *_5, 065 *0.105* SRI P0S*0.0*4JL) 

C^LL SYMBOL { 1 .'65t*4 .855,0 . 105 .SRlRAl *'b .'0 *37 ) 

IF( IMAX.EQ.I) GO TO 64020 

CALL SYMBOL(.95,4.6X.0.1225illHRE:FLECTOj^ 2*0,0*11) 

CALL SYMBOLd .05*4.40 , 0. 105*H2AMr’A, 0 ._0_*_18J 

CALL SYMDOL(1.05*4^19*0.105*SR2PoS*0.0*4i) 

CALL Symb ol ( i , o 5_*^3 . g 6_* o . i 05 *sr2 R aj ,_o_. o«37) 

IF( IMAX.E0.‘2) go To' 64020 

CALL 5YMH0L< .95*3, 735*0. 1225tliHREFLECT0R 3*0.0*11) 

CALL i‘YMBOL(1.05'.3.525*0.105*R3AMPA.0;oVl8) 

CALL SYMBOL n . 05 1 3 . 3 1 5* 0 , 1 05 * SR 3 POS * 0 . 0 :*i ) 

'call SyMBOLd. 05, 3. 105,0. 105, SK5RAT,0, 1)^37) 

64020 CALL SYMBOL( ,95,2.72, 0.14, 19HSCATTERED MULT I PATH* 0 , 0 * 19 ) 

CALL ~SYMSo'L(l.b5,i,51,6,105,SMA*0‘V6*l8') 

CALL SYMBOL) .95, 2, 125*0. 14*24HFR0NT-END RECEIVER NO ISE * 0 . 0 * 24 ) 

CALL SYMBOL) I ,05*1. '915, 0.105,3Hb“ =,0.0 t3f 

CALL PLOTd.l ,2.02,_3) 

CALL PLOf(l,15*2.0'2,2) ' 

CALL SYM 8 QH l..,14d.._?l_5*.P .07 *.1HN_*__0_. 0 *,1J 

CALL SYMBOLd.4* i. 915 * 0.105* SNA, 0 . 0 * 6 ) 

__CALL_S YMBOL d . 05 *J. . 705 , 0 . 1 0 5 * SNRGB A * 0 , 0 , l 6 )_ 

CALL SYMBOL (. 75 V 1 . 26570. 14, 20HCANDiOATE RECEIVER':. *0.0,20) 

CALL SYMBOLd. 05*1. 075 , 0*1D5*NKCVR*0, 0*50] _____ 

'cAL'l'sYMBOL (1 ,05*0 .67* 0.l05fNRCVR2,0.0 *30) 

CALL PLOT(D , 0 * - 0 41j “.3 ) 

GO TO)64100i642o6, 64300 ILBL 

_C 

c 

_1700P CALL ENDPLT 

STOP 

END 


APPENDIX B 


COMPUTER PLOTS 


AV. RMPLITUDE OF 
COMPOSITE SIGNRL 
RT TIME OF DIRECT 
SI&NRL CENTROID 
[SOLID LINE) 

RMS NOISE LEVEL 
tORSHED LINE) 


ESTIMRTION ERROR 
m/C RNCLE-ESTJ 
[DECREES) 


^RECEDING PAGE BLANK NOP KOtSa 



yK 


0.40 


0.80 1.20 l.GO 2.00 2.40 2.80 

TIME (SEC) SINCE STRRT QF FIRST SCAN 



/ V 


FIGURE B-l.a 

LQNC-TERM PLOT FOR 27 SCANS 
DATE: 10/07/7E 

PROCRRM: HLSRCVR 
JOBNAME: MLSRCMP 

S. H. IRWIN, JR. 


CHANNEL: A.ZIMUTH 

SIGNAL MAKEUP 

POSITION SPECIFIED 8'^ R^NEE IM NfiUTICflL MILES/ 

DEGREES OF AZIMUTH- DECREES OF ELEVATION, 

RATES SPECIFIED 8 VELOCI'Y IN KNOTS- DEGREES 
fiZIHUTH/SEC- DEGREES ELEVfiTION/SEC. 

DIRECT PATH 

AMPLITUDE = 1.00 

INITIRL R/C POSITION = 10.00, 30.00, 3. 

INITIRL R/C RRTE = -300.00, .00, .0 

SPECULAR MULTIPATH -NONE 


SCATTERED MULTIPATH 

AMPLITUDE = . 00 

FRGNT-END RECEIVER NOISE 

On= .10 
(SNR = 20.00 03) 

CANDIDATE RECEIVERS 

SQUARE GATE TRACKING RECEIVER 
(GO-BIT WORD LENGTH) 

ESTIMATION ERROR 

RMS ESTIMATION ERROR = .015 DEGREES 

RMS ERROR SPECIFICATION = .010 DEGREES 

^ SHORT-TEW PLOT OF COHPOSnE SIGNAL AVAILABLE 
O DETAILED SHORT-TERM PERFORMANCE PLOT HVf5lLHBLE 






scnN 
TK - 


'jcnN 

TK = 


scm 

TK = 


SCAN 
TK = 


SCAN 
TK = 


SCAN 
TK = 


1 

0 MS 


o 

o 

• 

o 


5 

300 MS 


3 

eoo MS 


13 

300 MS 


a 

o 

oj 


o 

o 

» ■ 

o 

CD 

CD 

OJ 


CD 

CD 

• ■ 

o 

o 

a 

OJ 


o 

o 

m • 

CD 


O 

o 

fsl' 

17 

1200 MS ^ 

o 

CD 


O 

o 

oj 

21 

1500 MS _ 





frj 





14.00 y.50 5.00 5.50 6.00 6.50 7.00 7.50 8.00 8.50 3.00 3.50 10.00 10. 5G 

RECEIVER SCfiN TIME (MILLISECONDS) 'l 




FIGURE B-Lb 

RMPLITUDE OF COMPOSITE SICNRL. PLUS NOISE 
DRTE: 10/07/75 

PROGRRM! MLSRCVR 
JQBNflME: MLSRCMP 

S, H. IRNIN. JR. 



CHANNEL: RZIMUTH 

SIGNRL MAKEUP 

posnim spEcn^iED b wwk im Msunca JtaEs. 

DECREES OF HZIHBTH- DECREES OF ELEVSTIUff, 

RfllES SFECIFJEa B ViUJOTTf IH tMOTS, DECREES 
(BIKUTif/SFC. AGREES ELEVOnOH/SEC. 

DIRECT PATH 

AMPLITUDE = l.OD 

INITIAL fl/C POSITION - 10.00, 30.00, 3,0D 

INITIAL n/C RATE = -300.00. .00. .0 

SPECULAR MULTIPATH -NONE 




SCATTERED MULTIPATH 

HriPl TTUDE - .00 

FRONT-END RECEIVER NOISE 

Oit= .10 

(SNR = 20.00 QB) 

CANDIDATE RECEIVER: 

SQUARE CATE TRACKING RECEIVER 
(GO-BIT WORD LENGTH) 



5CR!vl RNGLE 
(DEGREES) 

(ORSHED) 


CRNDIORTE 
RCVR OUTPUT 
(DEGREES) 

DIRECT PRTH 
COHPONENT ONLY 


RMPLITUDE OF 

MULTIPRTH 

INTERFERENCE 


RMPLITUDE OF 
COMPOSITE SIGNRL 
PLUS NOISE 


CRNDIDRTE 
RCVR OUTPUT 
(DEGREES) 

COMPOSITE SIGNAL 
PLUS NOISE 








y.oo 


U.50 


5,00 


5.50 G.OO 6.50 7.00 

RECEIVER SCAN TIME (MILLISECONDS) 


7.50 


8.00 


8.50 


8.00 


9.50 10.00 10.51 



FIGURE B- t.c 



o 


AMPLTTyOE OF 
Dlf^ECT PRTH 
COMPONENT 
[SOLID LINE) 


SHORT-TERM PLOT FOR SCAN 17 
DATE; 10/07/75 
PROGRAM; ML5RCVR 
JOBNRME; MLSRCMP 

S. H. IRNIN, JR. 


CHANNEL; AZIMUTH 


SIGNAL MAKEUP 

FOSinO!! SPEtiriED B ' RRXCE IN WUTItlS. BILE5. 

OEGHEES (F (UIHUTH. D‘'GR£ES OF ELEVRTMS). 
miES SPECIFIED Bt VELOcnV II! KltOIS. DECREES 
flZIKUlHFSEC, DEGREES ELEVBnOH/SEC 

( 30.00 DEG.) DIRECT PRTH 

RMPUTUDE = i.DO 

INITIRL fl/C P3SITI0N = 10.00. 30.00. 3.00 

INITIfiL R/C RATE = -300.00. .00. .0 

SPECULAR MULTIPATH -NONE 




SCATTERED MULTIPATH 

AMPLITUDE = .00 

FRONT-END RECEIVER NOISE 

On= .10 
(SNR = 20,00 OB) 

CANDIDATE RECEIVER: 

SQUARE CATE TRACKINC RECEIVER 
160-BIT WORD LENGTH] 






R/C RT 30.00 DEGREES AZIMUTH 
TK = 1200MSEC. 




RV. RMPLITUDE 
COMPOSITE SIGNAL 
AT TIME or DIRECT 
SIGNAL CENTROID 
(SOLID LINE) 

RMS NOISE LEVEL 
(DASHED LINE) 


-! . ... I I i f 


O 

CD 



TIME (SEC) SINCE START OF FIRST SCAN 


ESTIMATKDN ERROR 
(R/C ANGLE-EST.) 
(DEGREES) 




FIGURE B-2.Q 

LONG-TERM PLOT FOR 27 SCANS 
DATE: 10/08/75 

PROGRAM: MLSRCVR 
JOBNAME: MLSRCAZ 

S. H. IRNIN. JR. 

CHANNEL: AZIMUTH 

SIGNAL MAKEUP 

POSU'ON SPFCIPIED BY RANGE IN NflUTICfiL raLES. 

DEGREES OF fiZIHUTH^ DEGREES OF ELEVRTION. 

RATES SPB:IFIED BY VEtOCIT'' IN KNOTS, DEGREES 
AZIHUTH/SEC, DEGREES ELEVHTION/SEC. 

DIRECT PATH 

RHPLITUDE = 1.00 

INITIfiL R/C POSITION = 10.00, 30.00, 3.00 

INITIRL R/C RRTE = -300.00, .00, .0 

SPECULAR MULTIPATH -NONE 




1 

I 

! 


SCATTERED MULTIPATH 

RHPLITUDE = . 00 

FRONT-END RECEIVER NOISE 

On= .140 
tSNR = 8.00 DB) 


CANDIDATE RECEIVER: 

SQURRE GRTE TRflCid'NC RECEIVER 
(GO-BIT W(3R0 LENGTH) 


ESTIMATION ERROR 

RMS ESTIMATION ERROR = .053 DEGREES 

RHS ERROR SPECIFICRTION = .010 DEGREES 


5^ SHDRT-TERH PLOT OF CflHPOSITE SIGNAL AVAILABLE 
O DETAILED SHORT-TERM PERFORMANCE PLOT AVAILABLE 






•VY 



/ 








'fW^I 




nLjliii!.UK\.rnitlllfi l^LJu.ti.l Ikl^.ill J ai 










2.00 2.50 3.00 3.50 4,00 4,50 5.00 



5.00 


7.50 


8. 00 


8.50 


9.00 


5.50 G.OO S.50' 7.00 

RECEIVER SCAN TIME (MILLISECONDS) 


slso 


10.00 


10.50 11.00 113 


i 

) 

1 


gOUJOKgJ®AMB 
' ( 



! 


1 


i 


FIGURE B“2.b 

RMPLITUDE OF COMPOSITE 5IGNRL PLUS NOISE 
DATE: 10/08/75 

PROGRRM; MLSRCVR 
JOSNflME: MLSRCR2 

S. H, IRNIN, JR. 

CHRNNEL: RZIMUTH 

5ICN9L MRKEUP 

POSIIION SPFCIFIED BY lf(W!E IN IHUfTIML fUlfS, 

DECfiEES OF BZIMU1H. KGREE5 OF aEVRTIOJ). 

RSIES SFECIPIEQ 6' YaOCnr IN SHOTS. 0EGHEE5 
OIIHUTH/SEC, DEGFCE5 ELEmUON/SEC, 

DIRECT PRTH 

RHPLITUDE = 1.00 

IMITIRL R^C POSITION = 10.00, 30.00, 3.00 

INITIRL fl/C RATE = -300.00, .00, .0 

SPECULAR MULTIPATH -NONE 



SCATTERED MULTIPATH 

AMPLITUDE = .00 

FRONT-END RECEIVER NOISE 

On= .110 
(SNR = B.OO DB) 

CANDIDATE RECEIVER: 

5QURRE CflTE TRflCKINC RECEIVER 
(GO-BIT WORD LENGTH) 










3CRN RNGLE 
(DEGREES) 

(DRSHED) 


CRNDIDflTE 
RCVR OUTPUT 
(DEGREES) 

DIRECT PRTH 
COMPONENT ONLY 


RMPLITUDE OF 

MULTIPATH 

INTERFERENCE 


AMPLITUDE OF 
COMPOSITE SIGNAL 
PLUS NOISE 


CANDIDATE 
RCVR OUTPUT 
(DEGREES) 

COMPOSITE SIGNAL 
PLUS NOISE 






i bToo 5! 50 gToo gIso tToo tTso Soo aToo also io.ro 16.50 li.Go ] 

RECEIVER SCRN TIME (MTLLISECOND 33 





i 



rsi 

AMPLITUDE 0F 

FIGURE B-2.C 


DIRECT PfiTH 

SHORT-TERM PLOT F(3R SCAN 17 

a 

• 

COMPONENT 

DATE: 10/08/75 

o 

(SOLID Llf€) 

PROGRAM: MLSRCVR 


JQBNRME: MLSRCRZ 

3. H. IRWIN, JR. 


( 30.00 DEG.) 


CHRNNEL: AZIMUTH 

SIGNAL MfifCEUP 

POSnim 5PEOTJB1 BY RfflKE IM MaUIICfiL KIL£S- 
DEGREES OF flZIKUTH- DEGfiEES OF ELEVflTIOH. 

RflTES SPECIFKO B VELOCITY IN RtiSIS. DEGREEB 
BZIKUTH/SEC- DECREES ELEVOTION/SEC. 

DIRECT PATH 

RHPLTTUDE = i.OD 

INITIAL fl/C POSITION = 10. OG^ 30.00. 
INITIRL fl/C RATE = -300.00. .00. 


SPECULAR MULTIPATH -NONE 





{ 


23.30 DEG.) 


SCATTERED MULTIPATH 

AMPLITUDE = .00 

FRONT-END RECEIVER NOISE 

On= .40 
CSNR = 8.00 D0J 

CANDIDATE RECEIVER: 

SQUARE GATE TRACKING RECEIVER 
CGC-BIT WORD LENGTHS 


q 


10.00 10.50 


11.00 11.50 12.00 


12.50 


A/C AT 30.00 DEGREES AZIMUTH 
Ti; = 1200MSEC. 


AV. fiMPLITUDE OF 
COMPOSITE SIGNRL 
RT TIME OF DIRECT 
SICMRL CENTROID 
(SOLID LINE) 

RMS NOISE LEVEL 
(DRSHED LINE) 


.1 ...... - ..a. I 


o 

CD 



TIME (SEC) SINCE STRRT OF FIRST SCRN 



FIGURE B-3.a 

LONG-TERM PLOT FOR 27 SCRNS 
DRTEs 10/07/75 
PROGRAM; ML3RCVR 
JOBNRME; MLSRCJO 

S. H. IRNIN, JR. 

CHANNEL; AZIMUTH 
SIGNAL MAKEUP 

POSITION SPECIFIED B BRNGE IN NAUTICAL MILES, 

DEGREES OF AZIMUTH, DEGREES OF ELEVATION. 

RATES SPECIFIED BY VaOCHY IN KNOTS, DEGREES 
AZIHUTH/SEC, DEGREES ELEVATION/SEC. 

DIRECT PATH 

AMPLITUDE = 1.00 

INITIAL R/C POSITION = 10.00, 30.00, 

INITIAL R/C RATE = -300.00, .00, .0 

SPECULAR MULTIPATH 

REFLECTOR 1 

AMPLITUDE = .80 

INITIAL POSITION = 1.00, 33.00, 1.85 

INITIAL RATE = .00, -3.04, .00 

R.R PHASE DIFFERENCE = IS0“ AT PULSE 
COINCIDENCE 


SCATTERED MULTIPATH 

AMPLITUDE = .00 

FRONT-END RECEIVER NOISE 

On= .40 
(SNR = 8.00 D3) 

CANDIDATE RECEIVER; 

SQUARE CATE TRACKING RECEIVER 
(GO-BIT WORD LENGTH) 

ESTIMATION ERROR 

RMS ESTIMATION ERROR = .24,7 DEGREES 

RMS ERROR SPECIFICATION = .010 DEGRFES 

^ SHORT-TERM PLOT OF COMPOSITE SIGNAL AVAILABLE 
0 DETAILED SHORT-TERM PERFORMANCE PLOT AVAILABLE 









00 5.50 G.OO e.50 7.00 

RECEIVER SCRN TIME (MILLISECONDS] 


8.00 8.50 3.00 3.50 10,00 10.50 11.00 11.50 . 















FI&URE B-3.b 

AMPLITUDE QF COMPOSITE SIGMRL PLUS NOISE 
DATE: 10/07/75 

PROGRAM: MLSRCVR 
JGBNRME; MLSRCJO 

3. H. IRWIN, JR. 



CHANNEL: AZIMUTH 

SIGNAL MAKEUP 

ii35nim sittimo e Rfws w t«uiica miles. 
PEEEIS (F fiZDUnH. OIEfiHS (F ELECTION, 

IttTEB SPECIFIED B ' VaOClTV IM KWTS, DESEES 
(antUTH/SEC eiBHEES aEVHTKN/sa: 


DIRECT PATH 

AMPLITUDE = 1.00 

INITIRL fi/C POSITION = 10.00. 30.00. 3.C0 

INITIRL fl/C RATE = -300.00. .00. .0 

SPECULAR MULTIPATH 

RFFLECTOR 1 

RMPLITUDE = .80 

INITIRL POSnrClN = 1.00. 33.00. 1.85 

INITIAL RATE = .00. -3,04. .00 

RF. PHASE DIFFERENCE = 100“ At PULSE 





SCAN RNGLE 
[DECREES) 

CDfiSHED) 


CfINDIDRTE 
RCVR 0UTPUT 
(DEGREES) 

DIRECT PATH 
COMPONENT ONLY 


AMPLITUDE OF 

MULTIPATH 

INTERFERENCE 


AMPLITUDE OF 
COMPOSITE SIGNAL 
PLUS NOISE 


CANDIDATE 
RCVR OUTPUT 
(DEGREES) 

COMPOSITE SIGNAL 
PLUS NOISE 
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AMPLITUDE 0F 
DIRECT PATH 
COMPONENT 
(SOLID LINE) 


30.00 QEQ.) 


20.33 DEC.) 


iO.OO 10.50 11.00 11.50 12.00 12.50 


y --X'. 


Fl&URE B-3.C 

SHORT-TERM PLOT FOR SCRN 17 
DRTE: 10/07/75 

PR(3GRHM: MLSRCVR 
JOBNRME; MLSRCJO 

S. H. IRWIN, JR. 


CHRNNEL: RZIMUTH 

SIGNAL MAKEUP 

rasmofl SPECIFIED BT RfflK IN NHUIICffl. HO£5- 
0E£f)£E3 OF A2IHU1H. DEtKSS OF ELEVBIIIW. 

){{)i£S jsiEciFiED er vELocnr in knots, oegsees 
nZIHUTH/SEC. degrees EUVflTlON/SEC, 

DIRECT PATH 

AMPLITUDE = l.DD 

INITIAL A/C POSITION = 10.00. 30.00. 3.0D 

INITIAL A/C RATE = -300.00. .00. .0 

SPECULAR MULTIPATH 

REFLECT0R 1 

AMPLITUDE = .00 

INITIAL POSITION = 1.00. 33.00. 1.85 

INITIAL RATE = .00. -3. Oil. .00 

RF PHASE DIFFERENCE = ISO® AT PULSE 
COiNClDENCE 


SCATTERED MULTIPATH 

AMPLITUDE = .00 

FRONT-END RECEIVER N(3I3E 

On= .40 
(SNR = 8.00 DB) 

CANDIDATE RECEIVER: 

SQUARE GATE TRACKING RECEIVER 
(GO-BIT WORD LENGTH) 


fl/C AT 30.00 DEGREES AZIMUTH 
TK = 1200MSEC. 


I. 


RV. RMPLITUDE 0F 
COMPOSITE SIGNRL 
RT TIME OF DIRECT 
SIGNRL CENTROID 
(SOLID LINE) 

RMS NOISE LEVEL 
(DRSHED LINE) 



o 

o 


*^.00 


ol^o olso ITio ITgo Zoo 

TIME (SEC) SINCE STRRT OF FIRST SCRN 


CD 

rsj 

n 

CD 

00 

» 

o 


ESTIMRTION ERROR 
(R/C RN(GLE-EST.) 
(DE(;REES) 


o 




FIGURE B-H.a 

LONG-TERM PLOT FOR 27 SCRNS 
DRTE: lO/OG/75 

PROGRAMS ML3RCVR 
JGBNRME: ML3RCFG 

3, H. IRWIN, JR. 


CHRNNEL: AZIMUTH 

SIGNAL MAKEUP 

PQSniON SPECIFIED BY RflMSE IN NflUTICflL MILES, 

DEGREES OF flZIKUTH, DEGREES OF ELEVRTIDN. 

RRTE5 SPECIFIED BY VELOCnV IN KNOTS, DEGREES 
flZIHDTH/SEC, DEGREES ELEVRTI0N/SEC. 

DIRECT PATH 

AMPLITUDE = 1.00 

INITIAL R/C POSITION = 10.00, 30.00, 3.00 

INITIAL R/C RATE - -300.00, .00, .0 

?.80 


SPECULAR MULTIPATH 

REFLECTOR 1 
AMPLITUDE = .80 

INITIAL POSITION = i.OO, 33.00, 1.05 

INITIAL RATE = .00, -3.04, .00 

RF PHASE DIFFERENCE = 180° AT PULSE 
COINCIDENCE 


SCATTERED MULTIPATH 

AMPLITUDE = . 00 

FRONT-END RECEIVER NOISE 

Oh= .10 
(SNR ’ 20.00 DB) 

CANDIDATE RECEIVER^ 

SQUARE GATE TRACKING RECEIVER 
(GO-BIT WORD LENGTH) 

ESTIMATION ERROR 

RMS ESTIMATION ERROR = .203 DECREES 

RMS ERROR SPECIFICATION = .010 DEGREES 


^ SHORT-TERM PLOT OF COMPOSITE SIGNRL flVfilLflBLE 
0 DETfilLL’D SHORT-TERM PERFQRHflNCE PLOT fiVfllLHBLE 







i -i, • . : 



FIGURE B-^.b 

AMPLITUDE OF COMPOSITE SIGNAL PLUS NOISE 
DATE: lO/OG/75 

PROGRAM: MLSRCVR 
JQBNAME: MLSRCFG 

S. H. IRWIN. JR- 

CHANNEL: AZIMUTH 

SIGNAL MAKEUP 

P053T3OM SPECIFIED Br RTOE IH (CTJTlCfiU >m£Si 
DEGREES CF fiMHUIH- DEGREES IF ELEVSTOTI. 
fWTES SPECIFIED Br DEUJCm IN KHBTS. DECREES 
flZIHUTH/SEC- DECREES aEVSTTON/SEC, 

DIRECT PATH 

RMPLITUDE = 1.00 

INITIRL R/C POSITION r 10.00. 30.00. 3.00 

INITIAL n/C ARTE = -300.00- .00. .0 

SPECULAR MULTIPATH 

REFLECTOR 1 

AMPLITUDE = . BO 

INITIAL POSITION = l.OD, 33.00. 1.05 

INITIRL RATE = ,00, -3. OH. .00 

RF PHASE DIFFERENCE - 180“ AT PULSE 
COINCIDENCE 



SCATTERED MULTIPATH 

AMPLITUDE = .00 

FRONT-END RECEIVER NOISE 

Um= .10 
:SNR = 20,00 OB) 

CANDIDATE RECEIVER: 

SQUARE GATE TRACt^ING RECEIVER 
I GO-BIT WORD LENGTH) 


. , i 
1 



SCHN ANGLE 
(DEGREES) 

(DASHED) 


CANDIDATE 
RCVR OUTPUT 
(DEGREES) 

DIRECT PATH 
COHP0NENT ONLY 


AMPLITUDE OF 

MULTIPATH 

INTERFERENCE 


AMPLITUDE OF 
C(3MP0SITE SIGNAL 
PLUS NOISE 


CANDIDATE 
RCVR OUTPUT 
(DEGREES) 

COMPOSITE SIGNAL 
PLUS NOISE 





\ foi rom, FRAME 


} 


J 


... : • -'J-. 




a . t 


Fl&URE B-^.c 

SH0RT-TERM PLOT FOR SCRN 17 
DATE: iO/OG/75 

PROGRAM: MISRCVR 
JOBNAME: MLSRCFG 

S. H. IRWIN, JR. 



CHANNEL: AZIMUTH 

SIGNAL MAKEUP 

POSniDM SPECIFIED BV RRMCE IN NaLFTICfiL WILES. 

DECTEES OF flllWUTH. DESHEES OF ELEVflTIOH. 
fIflTES SPECIFIEO BY VELOCnv IN KNOTS. DECREES 
q2lHmH/£EC. DECREES EL£VflTION/SEC. 

DIRECT PATH 

AMPLITUDE = 1.00 

INITIAL B/C POSITION = IQ.OO, 30. DD, 3.00 
INITIAL A/C RATE = -300.00, .00, .0 

SPECULAR MULTIPATH 

REFLECTOR 1 

AMPLITUDE = .00 

INITIAL POSITION = l.DQ, 33.00, 1.05 

INITIAL RATE = .00, -3.011, .00 

RF PHA5E DIFFERENCE = 180° AT PUU3E 
COINCIDENCE 



SCATTERED MULTIPATH 

AMPLITUDE = .00 

FRONT-END RECEIVER NOISE 

Oh= .10 
(SNR = 20.00 OB) 

CANDIDATE RECEIVER: 

SQUARE GATE TRACKING RECEIVER 
(SO-BIT WORD LENGTH) 

A/C AT 30.00 DEGREES AZIMUTH 
TK = 1200MSEC. 



RV RHPLITUDE OF 
COMPOSITE SIGNRL 
RT TIME OF DIRECT 
SIGNRL CENTROID 
(SOLID LINE] 

RMS NOISE LEVEL 
(DRSHED LINE) 



ESTIMRTION ERROR 
(R/C RNGLE-EST .. ) 
(DECREES) 







FIGURE B - 5.a 

LONG-TERM PLOT FOR 21 SCPNS 
DRTEs 10/10/75 
PROGRRM' MLSRCVR 
JOBNRMEs HLSRC5C 

S. H. IRNIN, JR- 


CHRNNEL: RZIMUTH 

SIGNRL MRKEUP 

POSniOK SPECIFIED BY RfSNGE IN NflUTICRL MILES 
DEGREES OF flZIHUTH/ DEGREES OF ELEVATION. 

RATES SPECIFIED BY VELOCnY IN KNOTS. DEGREES 

nziHUTH/sa:. degrees elevation/sfc. 

DIRECT PRTH 

fiMPLITUDE =: 1.00 

INITIRL R/C POSITION = 10.00, 30.00, 3.00 

INITIAL R/C RRTE = -300.00. .00. .0 

SPECULAR MULTIPP.TH 

REFLECTOR 1 

RNPLITUDE = . BO 

INITIRL POSITION = 1,00 33.00 1.85 

INITIRL RRTE = .00 -3.011, .00 

RF PHR5E DIFFERENCE = o'* RT PULSE 
COINCIDENCE 


SCATTERED MULTIPATH 

AMPLITUDE = .00 

FRONT-END RECEIVER NOISE 

CTn= .10 
(SNR = 20.00 DB) 

CANDIDATE RECEIVERS 

SQURRE GATE TRACKING RECEIVER 
(GO-BIT WORD LENGTH) 

ESTIMATION ERROR 

RMS ESTIMATION ERROR = ,1G2 DEGREES 

RMS ERROR SPECIFICATION = .010 DEGREES 

^ SHORT-TERM PLOT OF CWOSHE SIGNAL flVRILABLE 
O DETAILED SHORT-TERH PERFORMANCE PLOT fiVATLASLE. 




L.OO 
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FIGURE B-5.b 

RMPLITUDE OF COMPOSITE SIGNAL PlUS NOISE 
DATE: iO/iO/75 

PROGRAM; MLSRCVR 
JGBNAME: MLSRC5C 

S. H. IRWIN, JR. 

CHANNEL; AZIMUTH 
SIGNAL MAKEUP 

FBsmm SFECtf IH] BT HfW!£ IM HSUTTCftL MJL£S- 
DEG«E£S OF BHHinil, DEGREES OF ELEVRTTCN. 

RATES SPECIFIS3 BY VaoCHY IH WIOTS. BEGREES 
RZIKlfm/SEC< DECREES ELEVRTimt^C. 

DIRECT PATH 

AMPLITUDE = l.OD 

INITIAL fl/C POSITION = 10. DO, 30.00, 3.00 

INITIAL fl/C RATE = -300 00, .00, .0 

SPECULAR MULTIPATH 
REFLECTOR i 
RMPLITUDE = .00 

INITIAL POSITION = 1,00, 33.00, 1.B5 

INITIAL RATE = .00 -3.0A .00 

RF PHASE DIFFERENCE = O" RT PULSE 
COINCIDENCE 



SCATTERED MULTIPATH 

A'tPLlTUDE = .00 

FR0NT-END RECEIVER NOISE 
On= .10 
(SNR = 20.00 OB) 

CANDIDATE RECEIVER; 

SQUARE GATE TRACKING RECEIVER 
(GO-BIT WORD length; 





SCRN RNCLE 
(DEGREES) 

(DR3HED) 


CRMDIDRTE 
RCVR OUTPUT 
(DEGREES) 

DIRECT PRTH 
C(DMP0NENT ONLY 


RMPLITUDE OF 

MULTIPATH 

INTEREERENCE 



AMPLITUDE OF 
COMPOSITE SIGNAL 
PLUS NOISE 


CANDIDATE 
RCVR OUTPUT 
(DEGREES) 

COMP(33ITE SIGNAL 
PLUS NOISE 


o 



0.00 0.50 1.00 1.50 2.00 2.50 3.00 3.50 4.00 4.50 5.00 








00 5.50 G.OO G.50 7.00 7.50 8.00 8.50 9.00 9.50 10.00 10.50 11.00 11.50 

RECEIVER SCON TIME (MTLLTSEC0NOS) 


gx)ijx)ua; gRAKs ^ 

V 






CD 

O 
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o 


AMPLITUDE OF 
DIRECT PATH 
COMPONENT 
CSOLID LINE) 


FIGURE B-5,c 

SHORT-TERM PLOT FOR SCAN 17 
DATE: 10/10/75 

PROGRAM: MLSRCVR 
JOBNAMEs MLSRC5C 

S. H. IRWIN, JR. 



[ 


30.00 DEG.) 


CHANNEL: AZIMUTH 

SIGNAL MAKEUP 

POSnlBN SPECIFIED BV W'lQE ID NflUlTCRL M’LES 
DECREES OF R21MUIM, DECREES OF ECtVRItBN. 
times SPECIFIED B-i VELDcm IN KIHIS. DH31EES 
flllKnH/SEC. decrees ELEVmim/SEC. 

DIRECT PATH 

fiMPLlTUDE - 1.00 

INITIAL R/C POSITION = 10-00, 30.00, 

INITIAL A/C RATE = -300.00 .00 

SPECULAR MULTIPATH 

REFLECTOR 1 

AMPLITUDE - .eO 

INITIAL POSITION = 1.00 33.00 1 

INITIAL RATE ' .00, -3.0H, .00 

RF PHASE DIFFERENCE = O' AT PULSE 
COINCIDENCE 




V 


I 


29.86 DEG.) 


SCATTERED MULTIPATH 

AMPLITUDE = .00 

FRONV-END RECEIVER NOISE 

CTn= .10 
tSNR = 20.00 OBI 

CANDIDATE RECEIVER: 

SQUARE CATE TRACKING REC.”VER 
lEO-BlT WORD LENGTH) 


^ iLoo li.50 laToo ?2 . 50 


3.00 

.0 


.85 


10.00 10.50 


A/C AT 30.00 DEGREES AZIMUTH 
TK = 120GMSEC. 


1 


J 




i ... -..^-.-1 i 


RV. RMPLITUDE OF 
COMPOSITE SIGNAL 
AT TIME OF DIRECT 
SIGNAL CENTROID 
(SOLID LINE! 

RMS NOISE LEVEL 
(DASHED LINE) 


ESTIMATION ERROR 
(R/C RNGLE-ESTJ 
(DEGREES) 



o 

(JD 








FIGURE B-6.a 

LONG-TERM PLOT FOR 27 SCRNS 
DRTEs 10/10/75 
PROGRAM! MLSRCVR 
JOBNflME: MLSRC8E 

S. H. IRWIN, JR. 


CHANNEL! AZIMUTH 
SIGNAL MAKEUP 

PO'JITIOfJ SPECIFIED BY RANGE IN NflUTICflL MILES/ 

DEGREES OF AZIMUTH/ DEGREES OF ELEVATION. 

RATES SPECIFIED BY VELOCITY IN KHOTS/ DEGREES 
AZIHUTH/SEC/ DEGREES ELEVATIOM/SEC- 

DIRECT PRTH 

RMPLITUDE = 1.00 

INITIRL n/C POSITION = 10.00/ 30,00/ 3.00 

INITIRL R/C RRTE = -300.00/ .00/ .0 

SPECULAR MULTIPATH 

2.80 REFLECTOR 1 

RMPLITUDE = ,80 

INITIRL POSITION = 1.00, 33.00, 1.85 

INITIRL RRTE = .00/ -3.011/ .00 

RF PHASE DIFFERENCE = 90 ° RT PULSE 
COINCIDENCE 


2.80 


SCATTERED MULTIPATH 

AMPLITUDE = .00 

FRONT-END RECEIVER NOISE 

On = .10 

(SNR = 20.00 DB3 

CANDIDATE RECEIVERS 

SQUARE GATE TRACKING RECEIVER 
(GO-BIT WORD LENGTH) 

ESTIMATION ERROR 

RMS ESTIMATION ERROR = .173 DEGREES 

RMS ERROR SPECIFICRTION = .010 DEGREES 

^ SHQRT-TERK PLOT CF CewT'CSITE SIGJ'tAL AVAILABLE 
0 DETAILED SHORT-TERM PERFBRMAMCE PLOT AVAILABLE 








7.50 8. DO 8.50 9.00 3.50 
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F^6l^£ B-6i2 

£^_TT].trr OF !IS»‘§^rfE SISsSiL ^lE 
iHTE: 10/10/1S 

HL^C^ 

S. H. JR. 


CfiPJaCLi mMOlH 
SI^RL miBP 

Sf ««E 3 fitKESt SSS- 
3QCSF ijt3iiyiysffi«t:» 

»EjiCnrj6ia*^ UCTES 
ssssenu^:^ ££K<stBK.'S£. 

DIRECT pmn 

a»e>LITi£H = l.iS 

INITIJS- PESinSN = iO.EJ- 30.OJ- 3.ao 
ISITiaL fi/c RJ?TE = -3D3.3G, .00. .0 

5PECUL53R ?!ULTIP9TH 

RSTICira I 
RHPLITUDr = ,BD 

INITIflL PBSmON = 1.30. 33.03- 1.85 

INITIflL RHTE = .DD- -a.OU- .£K3 
RF FHB^ DIFFERENCE =- TO* «T PULSE 
COINOOENCE 








SCfiTTERED MULTIPATH 
fWPLrruoE = .00 

FRONT-END RECEIVER NOISE 

0(1= .ID 
(SNR = 20.00 db; 


CANDIDATE RECEIVER: 

SaaflRE CRTE TRRCiaNe RECEIVER 
(60-BIT WORD LENCTH) 




SCRN ANGLE 
[PEGREES] 

(DASHED) 


CANDIDATE 
RCVR 0UTPUT 
(DEGREES) 

DIRECT PATH 
C0MP0NENT ONLY 


AMPLITUDE 0F 

MULTIPATH 

INTERFERENCE 


AMPLITUDE OF 
COMPOSITE SIGNAL 
PLUS NOISE 


CANDIDATE 
RCVR OUTPUT 
(DEGREES) 

COMPOSITE SIGNAL 
PLUS NOISE 










AMPLITUDE OF 
DIRECT PRTM 
COMPONENT 
(SOLID LINE) 


FISURE B-6.G 

SHORT-TERM PLOT FOR SCAN 17 
DATE: 10/10/75 

PR0(JRAM: MLSRCVR 
JOSNAME: MLSRC8E 

S. H. IRWIN, JR. 



CHANNEL; AZIMUTH 
SICNRL MAKEUP 

PQSniCN SPEOFIEO BY RDHCE JM HHUTIGH, NILES. 

0EBNEE3 8F fflWnH. tlE«£S BF OEVflTaN. 

(WTES SPECIE BY IffiLOmY IN HNUTS. DEOffiES 
niNJIH/SEC. 0ECREE9 £UVflTI<M/X[:. 

DIRECT PATH 

RMPtlTUDE = i.DO 

INITIAL fl/C POSniON = iO,DO> 30.00, 3.00 

INITIAL R/C RRTE = -300.00, .00, .0 

SPECULAR MULTIPATH 

REFLEGT0R 1 

AMPLITUDE n ,80 

INITIAL POSITION = 1.00, 33.00, 1.65 

INHIRL RATE = .00, -3.0^, .00 

RF PHASE DIFFERENCE = AO* AT PULSE 
COINCIDENCE 



SCATTERED MULTIPATH 

RMPLITUDE = .00 

FRONT-END RECEIVER NOISE 

(Tn= .10 
(SNR = 20.00 db; 

CANDIDATE RECEIVER; 

SQURRE ERTE TRRCKlNE RECEIVER 
CGO^BIT WORD LENETH) 

A/C AT 30.00 DEf^REES AZIMUTH 
TK - 1200MSEC. 


a 





FIGURE B-7.Q 

L0Ne-TERM PLOT FOR 27 SCRNS 
DfiTEs 10/10/75 
PROORRMs MLSRCVR 
J0BNRME5 MLSRC7Q 

S, H. IRWIN, JR. 


CHRNNELs RZIMUTH 


SIGNRL NRKEUP 

P3Srtl3N SPECIFIED SV RfW6£ IN NflUTicaU MILES, 
DE6MEE3 0F RZIHUTH, DEEfiEES OF ELEVATION. 

RATES SPECIFIED 9V VELOCITV IN (tNOTS. DEGREES 
AZIMUTH/SEC> DECREES ELEVATION/SEG. 

DIRECT PRTH 

.WUTUDE = 1.00 

IMITIRL n/Z POSITION 10.00, 

INITIAL A/C RATE = ^300. 00, 

SPECULRR MULTIPRTH 


30.00, 

. 00 , 


REFLECTOR 1 

AMPLITUDE = .80 

INITIAL POSITION - 1.00, 33.00, 1.65 

INITIAL RATE = .00, -3.0U, .00 

RF PHASE DIFFIRENCE - 270* AT PULSE 
COINCIDENCE 


SCRTTERED MULTIPRTH 

AMPLITUDE = .00 

FRONT-END RECEIVER NOISE 

On= .10 
(SNR = 20.00 08) 

CRNDIDRTE RECEIVER; 

SQUARE CATE TRAGftINC RECEIVER 
leO-BIT WORD LENGTH) 

ESTIMRTI0N ERROR 

RMS ESTIMATION ERROR = .201 DEGREES 

RMS ERROR SPECIFICATION ~ .010 DECREES 

^ SHQHT-TERrt PLOT OF COMPOSTTE 5ICNHL flVfilLflBLE 
O DETAILED SHORT-TERN PERFORHHKCE PLOT AVAILAfiLE 


flV, (AMPLITUDE OF 
COMPOSITE SIGNAL 
AT TIME OF DIRECT 
SIGNAL CENTROID 
(SOLID LINE] 

RMS NOISE LEVEL 
(DASHED LINE) 



TIME (SEC) SINCE START OF FIRST SCAN 




ESTIMATION ERROR 
(R/C ANGLE-ESTJ 
(DEGREES) 








1 


wl^Sa|a*?tprn«^B 


SCflN 1 
TK = 0 MS 


SCAN i 
TK = 300 MS 


SCAN 3 
TK r GOO MS 


SCAN 13 
TK = 300 MS 


SCAN 17 
TK = 1200 MS 


SCAN 21 
TK - ISOO MS 







8.00 


8 


1:00 5.50 6.00 6.50 7.00 

RECEIVER SC8N TIME! MILLISECONDS 1 

t 

r 



ISAME*! 




7.50 
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{■ 

y 

? 





„JlA_ 




figure B- 7.b 

flHPLlTUDE 0F COMPOSITE SIGNAL PLUS NOISE 
DATE: 10/10/75 
PROGRAM: MLSRCVR 
JOBNAME: MLSRC7Q 

S, H, IRWIN. JR, 


CHANNEL: AZIMUTH 

SIGNAL MAILEUP 

pdsmw specifiQ) sr me m nautim. mms< 

OECHEEa OF nznMTtt/ aecii£E!t OF £L£vnT 10 H. 

(HIES sfwiFieo BV veuflcnF ih UNers. oecbees 
(CMSTH/seE. oemees euevBTjoM/sec. 

DIRECT PATH 

AMPLITUDE = 1.00 

INITIAL fi/C POSITION = 10.00, 30,00, 3.00 

INITIAL A/C RATE = -300.00, .00, .0 

SPECULAR MULTIPATH 
REFLECTOR 1 
AMPLITUDE = ,00 

INITIAL POSITION = i.GO, 33.00, 1.05 

INITIAL RATE = .00, -3. OH, .00 

RF PHASE DIFFERENCE = £T0® AT PULSE 
COINCIDENCE 







12.50 


SCATTERED MULTIPATH 

AMPLITUDE = .00 

FRONT-END RECEIVER NOISE 
Gr«= .10 
tSNR = 20.00 DB) 

CANDIDATE RECEIVER: 

SQUARE EHTE TRAdUNO RECEIVER 
teO-BIT WORD LENGTH) 









DIRECT PRTH 
COMPONENT ONLY 



AMPLITUDE OF 

MULTIPRTH 

INTERFERENCE 


RMPLITUDE OF 
COMPOSITE SICNRL 
PLUS NOISE 


CflNDIDflTE 
RCVR OUTPUT 
(DEGREES) 

COMPOSITE SIGNRL 
PLUS NOISE 





7.50 


8.00 


8.50 


9.00 


9.50 


5.00 5.50 G.OO G.50 7.00 

RECEIVER SeON TIME (MILLISECONDS) 


10.00 10.50 11.00 11.’ 










RMPLITUDE OF 
DIRECT PRTH 
COMPONENT 
[SOLID LINE] 



FIGURE B-7.C 

SHORT-TERM PLOT FOR SCRN 17 
DflTEt 10/10/75 
PROGRAM! MLSRCVR 
J0BNAME! MLSRC7Q 

S. H. IRWIN. JR. 


CHANNEL! AZIMUTH 
SIGNAL MAItEUP | 

pomieN ipnirip vt tsut ot mutich. nilcs, 
vata 9 iciiuni. er omtiw. 

ffnes sKcinZD bv vascny m (nois, otsnets 
flznunv^. oesma aBnrim«c. 

DIRECT PATH 

finPLITUDE = 1.00 

IMITIflL «/C POSITION = 10.00. 30.00- 3.00 

INITIAL A/C RATE = -300.00, .00, ,0 

SPECULAR MULTIPATH 

REFLECTOR 
AMPLITUDE = .00 

INITIAL POSITION = 1.00, 33.00, 1.85 

INITIAL RATE = .00, -3.0H., .00 

RF PHASE DIFFERENCE = 270“ AT PULSE 
CGINCinENCE 



SCATTERED MULTIPATH I 

AMPLITUDE = .00 

FRONT-END RECEIVER NOISE 

Dh= . 10 
tSNR = 20.00 00) 

CANDIDATE RECEIVER: 

SQUARE GATE TRACKING RECEIVER 
160-811 WORD LENGTH) 

A/C At 30.00 DECREES AZIMUTH 
TIC. = 1200MSEC. 



